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Abstract. The distribution function for the first eigenvalue spacing in the 
Laguerre unitary ensemble of finite rank random matrices is found in terms of 
a Painleve V system, and the solution of its associated linear isomonodromic 
system. In particular it is characterised by the polynomial solutions to the 
isomonodromic equations which are also orthogonal with respect to a defor- 
mation of the Laguerre weight. In the scaling to the hard edge regime we 
find an analogous situation where a certain Painleve III' system and its asso- 
ciated linear isomonodromic system characterise the scaled distribution. We 
undertake extensive analytical studies of this system and use this knowledge 
to accurately compute the distribution and its moments for various values of 
the parameter a. In particular choosing a = ±1/2 allows the first eigenvalue 
spacing distribution for random real orthogonal matrices to be computed. 



1. Introduction 

The Laguerre unitary ensemble (LUE„. a ) of random matrices is specified by the 
eigenvalue probability density function (p.d.f.) 

(1.1) P (X u ...,X n ) 

1 n 

Il e ~ A3A ? (A* -A,) 2 , Ai,...,A n e [0,00), 



\r 

3=1 l<j<k<n 



where 

^ ra 

(1.2) cm, n ■= —j IT r ( n " 3 + l)r(m - j + 2). 



mi 

3=1 



The naming relates to the fact that (jl.ip is the eigenvalue p.d.f. of complex Hermit- 
ian matrices X with measure invariant under unitary conjugation X i— > UXU~~ l , 
proportional to the generalised Laguerre form 

(1.3) (detX) a e- Tr(x) . 
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In multivariate statistics (jl . 1|) is realised as the eigenvalue p.d.f. for the complex 
case of the so-called Wishart matrices X = Y^Y. Here Y is an N X n (N > n) 
rectangular matrix of i.i.d. entries with distribution N[0, 1] -HN[0, 1]. In this setting 
a = N — n, and so a is naturally a non-negative integer. The spectrum of complex 
Wishart matrices has found recent application in studies of wireless communication 
systems [34], where the matrix Y consists of the complex amplitudes of various 
channels of transmitted waves as received by the antennas. 

The matrix structure Y^Y is relevant to the study of the eigenvalues of the 
(n + N) x (n + N) Hermitian matrix 

(1.4) X := 

Thus one has that X has in general N — n zero eigenvalues, with the remaining 
eigenvalues given by ± the positive square roots of the eigenvalues of Y'Y (see e.g. 
[5]). This matrix structure has application to the study of the Dirac equation in 
the context of quantum chromodynamics |36j . There most interest is in the scaling 
behaviour of the smallest eigenvalues. 

In the study of matrix Lie algebras one encounters antisymmetric matrices 
(X T = —X) with pure imaginary complex elements. Specifically, such matrices 
are the Hermitian part of the matrix Lie algebra 

(1.5) i x (so(n, C)) := {i times n x n skew symmetric complex matrices}. 

If the independent imaginary complex elements are i.i.d with distribution iN[0, 1], 
then the p.d.f. of the positive eigenvalues is proportional to 

m 

(1.6) n = 2m even: Jexp(-Af) {\\ — A|) 2 , 

j=l \<j<k<m 
rn 

(1.7) n = 2m + 1 odd: J \) exp(-Af ) (\\ - \]f . 

j — l l<j<k<m 

This ensemble will be denoted by AS(n). Under the change of variables A| i— > Aj 
these reduce to the LUE nitl with parameters a = — 1/2, 1/2 respectively. 

Antisymmetric Hermitian matrices X can be used to parameterise real orthog- 
onal matrices R with determinant +1 (and thus, by definition, members of the 
classical group + (n)). Thus we can write R according to a Cayley transformation 

Jl„ - iX 

Note from this that the property that the eigenvalues of X come in ± pairs is 
consistent with the property that the eigenvalues of R come in complex conjugate 
pairs e ±t6 . This can be used (see e.g. O) to show that with the matrix R chosen 
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with uniform (Haar) measure, the eigenvalue p.d.f for the eigenvalues with angles 
< < 7r is proportional to 

(1.9) n = 2m even: (cos Ok — cos Oj) 2 , 

l<j<k<m 
m 

(1.10) n = 2m +1 odd: ^[(l-cosfl,) (cos(9 fc - cos^) 2 . 

j—l l<j<k<m 

Note that for 8[ — > these have the same leading behaviour as (|1.6[) . (|1.7p with 
A; — * 0. This is consistent with the fact that the to — > oo scaled joint distribution 
function for the k smallest eigenvalues, p(u\ say, is the same for both ensembles, 

(1-11) 

m-oo^v/m/ ^ fc > V 2^' ' 2^/to^ m^ooVm/ ^( fe ) V TO ' TO 

where n = 2to, 2to + 1. This scaling is such that the average spacing between 
eigenvalues approaches unity as k — > oo. The distribution (jl.lljl is of primary 
importance in the study of the spectral interpretation of L- functions [5D] , [3D] . From 
the remark below (11.71) we know that 



(L12) (27^) • • • (^k) ■ ■ ■ - ^ = p^) Em ' a (^ ■ ■ • 

where for n = 2m, a = —1/2, while for n = 2m + 1, a = 1/2. Consequently 
1.13) hm — p' '{ 

m— >oo \TO/ W TO TO 

= i im f^Vxi x fcP LUE '-(^i ^) 

m— >oo \2m J 4to ' ' 4m 

Thus knowledge of the distribution P/^ Em '" for a = ±1/2 suffices to compute the 
scaled limit of p? k , ^ . 

In the case k = 1 a number of different characterisations of p^ Em,a , which is the 
distribution of the smallest eigenvalue in LUE mi<1 , are known. First, with n 1— ► n + 1 
in (|l.ip for convenience, the p.d.f. of the smallest eigenvalue is given by fixing one 
of the coordinates at x\, and integrating the remaining over [xi,oo). Thus 

(1.14) p^ E "^( Xl ) = - 1 e-^xl 



/OO p OO ^ 

dXi . . . / dA„ e-^A^Aj - x-!) 2 (A fe - A 

-1 ^1 7 = 1 l<7<fe<n 



J 2 - 



One has that 



(1-15) pS E - +1 -(*0 = -^ UB - +1 -(x 1 ), 
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where 

/CO POO 
dXf- J dA n+ u>(Ai ) ... ) A n+ i), 

is the probability (gap probability) that no eigenvalues are in the interval (0,t). 

For integer values of the parameter a in (jl.ip . £; LUE ",a (£) W as studied by orthog- 
onal polynomial techniques in [llj . where it was evaluated as an a x a determinant, 
and by the method of Jack polynomials in [10j . giving an a-dimcnsional integral 
form. In [32] (see also p3]), for general Re (a) > —1, it was expressed in terms of a 
fifth Painlcvc transcendent. Explicitly, it was found that 

(1.17) E WE ^(t)=ex P n*^-U v (s 
where Uy(s) satisfies the Jumbo-Miwa-Okamoto er-form of the Painlcvc V equation 

(1.18) {ta") 2 -(a- ta 1 + 2{a') 2 + (i/ + vi + v 2 + v^a'f 

+ 4(l/ + cr')(^i + a')(v 2 + cr')(^3 + cr') = 0, 

with parameters 

(1.19) isq = v\ = 0, v 2 = n + a, v% = n. 
Alternatively the conventional Painleve V parameters are 

(1.20) a=^a 2 , (5 = 0, 7 = -2n - a - 1, 8 = ~, 
and in terms of the Okamoto parameters they are 

(1.21) v 2 — vi — 0, V3 — vi = n + a, V4 — v\ — n, v% — V4 = a. 

Because the eigenvalue density is strictly zero for A < 0, the neighbourhood of the 
smallest eigenvalue is referred to as the hard edge, and is denoted by HE a . As 
is consistent with (|1.13p . a well defined limit of (|1.17p is obtained by the scaling 
1 1— > t/An and n — > 00. Thus |33j 

(1.22) lim E LVE "-"(t/4n) := £ HE «(i) = exp ( [ — U ln ,(s)) , 

n->oo \J Q S J 

where Uuy (t) satisfies the Jimbo-Miwa-Okamoto cr-form of the Painleve III' equation 

(1.23) (ta") 2 - Vl v 2 (a') 2 + a' (Aa 1 - l)(er- ta') - ^(vi - v 2 ) 2 = 0, 
with parameters 

(1.24) v x = v 2 = a, 
and subject to the boundary condition 
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Our interest in this paper is the distribution p^ Em,a and its hard edge scaled 
limit. Analogous to (|1.14p , we see from that 



(1.26) p 



LUE„ +2 ,, 

( 2 ) n\c n+2 ,n+2+a 



'(x 1 -x 2 ) 2 (x 1 x 2 ) a 



. dX n ]J e-^X^Xi - Xl f{Xi - x 2 f J] ( A fc - A i) 2 ' 

where x\ denotes the smallest eigenvalue and x 2 the second smallest eigenvalue. In 
[T2] . for a integer, this was expressed as an (a + 2) x (a + 2) determinant. In the 
hard edge scaled limit this gave 



1 LUE„ +2i0 ,Il X 2 HE 



4 X 2 



■i'l 



-x 2 /4 



x det 



x 2 - a?i 

, ^2 



[lj+2-k{y/X2j\ j=l,...,a 
k=l,...a+2 

(fe-i)/2 



I j+2 -k(\/x 2 - Xi) 



3=1,2 
fc=l,...a+2 . 



We seek a Painleve type characterisation of p.26|) and its scaled limit, valid for 
general Re (a) > — 1. 

One use of knowledge of Pt 2 -^ m ' a is the computation of the distribution of the 
spacing between the smallest and the second smallest eigenvalues. Denoting this 
distribution by A n ^ a for LUE„ + 2. a , we have 



(1.28) 



dx 1 p^ E " +2 ' a (x 1 ,x 1 +y), y G K + . 



Important to our subsequent working is a rewrite of p.26[) and (|1.28|) in terms of 
an integral of the form 

(1.29) D n { Xl ,x 2 )[w{X)] 

:=i fdX,... fdX n f[w(Xi)f[{\i - xx^X, - x a ) J] (A fc -A,) 2 , 

' J 1 Jl 1=1 1=1 l<j<k<n 

where / denotes the support of the weight w(X). We have 



(1.30) p (2) (x,x + y) 
1 



Cn+2,n+2+a 

and 
(1.31) 

A n Av) 



{n+ i)(x + y)- Xy 2 [x{x + y)] a Dnhyi _y)[A 2 (A + x + y) a e- x X>0 ], 



y 2 e» 



C n +2,n+2+a 



dt t a (t - yfe-^^D^-y, -y)[X 2 {X + t) a e- x X>0 } 
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These equations exhibit the occurrence of a deformation of the Laguerre weight 
(1.32) w(x;t) := x 2 (x + t) a e~ x , x G K+. 

This deformed weight actually interpolates between two Laguerre weights - when 
t — ► then we have the general parameter a + 2 case, whilst if t — ► oo in the sector 
— 7r < arg(t) < 7r we have the special parameter situation with an exponent of 2. In 
fact virtually all of our analysis can be carried over to the more general situation 
where the exponent 2 is an arbitrary complex parameter suitably restricted. 

We begin in Section 2 by revising appropriate results from orthogonal poly- 
nomial system theory and apply this to the particular deformed Laguerre weight 
(|1.32p . This allows us, in Section 3, to characterise the distributions p.30p and 
()1.31|) by a solution of the fifth Painlcve equation and its associated linear isomon- 
odromic system (see Proposition 3.2). Section 4 is devoted to the determinant 
evaluations of those distributions for positive integer values of the parameter a. We 
proceed in Section 5 to the study of the hard edge limits 

ft qq\ HE a I s. ,. 1 LUE„, ,J;i X 2 , 

( ) P(2) (^):= n lim ^P {2) 

(1.34) A a (z) := lim L A ^ a {±-). 

It is found that these scaled distributions can be characterised by the solution of a 
certain Painlcve III' equation and its associated linear isomonodromic system (see 
Propositions 5.2, 5.5 and Remark 5.3). In Section 6 this characterisation is used to 
obtain the high precision numerical values of statistical characteristics of A a (z) for 
various integer values of a and for the values a = ±1/2, the latter being relevant 
to (|1.6[) . (|1.7p with the change of variable A| i— ► Xj. Let p^Jxx, x%) denote the 
scaled distribution of the eigenvalues e lBl , e 102 [9\, 6*2 > 0) closest to the origin in 
+ (2n + 1) and + (2n) respectively With the scaling chosen so that the bulk 
density is unity, it follows from (| 1 . 13[) and (|1.33[) that 

(1.35) pf 2) (xi,x 2 ) = 4 7 r 2 a; 1 a;2P™ ±1/2 (7r 2 x2,7r 2 a;2). 
Consequently 



oo 



A (y) := / dxp (2) (x,x + y) 
Jo 

/>oo 

(1.36) =4tt 2 J dxx(x + y)p^ ±1/2 (TT 2 x 2 7 TT 2 {x + y) 2 ). 

We use our results for P™ ±1/2 to provide the high precision numerical values of 
statistical characteristics of A ± (y). 
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2. Orthogonal Polynomial System 



2.1. Semi-classical Orthogonal Polynomials. Consider the general orthogonal 
polynomial system {p n {x)}^ =li defined by the orthogonality relations 



(2.1) 



dx w(x)p n (x)x r 



< m < n 
h„ m = n 



with / denoting the support of the weight w(x). We give special notation for the 
coefficients of x n and x n_1 in p n (x), 



(2.2) 



Pn{x) = J n X n + 7„,lX™ 1 + 



The corresponding monic polynomials are then 



(2.3) 

It follows from (|2~T|) that 



K n (x) = —p n (x). 
In 



dx w(x)(p n (x)) 2 = 7„/i„, 



and thus for p n {x) to be normalised as well as orthogonal we set j n h n = 1. A 
consequence of the orthogonality relation is the three term recurrence relation 



(2.4) 



a n+1 p n+ i(x) = {x- b n )p n {x) - a n p n -i(x), n>l, 



and we consider the set of orthogonal polynomials with initial values p_i = 
and po = 7o- The three term recurrence coefficients are related to the polynomial 
coefficients by [31], [14] 



(2.5) 

along with 
(2.6) 



7«-i 

7n 



, 7n,l 7n+l,l . -. 

On = , n > 1, 



In 



ln+1 



bo = — ^i, ao = 0, 70,1 = 0. 
7i 



A well known consequence of (|2.4p is the Christoffcl-Darboux summation 



(2.7) 



^2pj(x)Pj(y) 

3=0 



\p n {x)p n - 1 {y) - p n -i(x)p n (y)] 



x~y 

Central objects in our probabilistic model are the Hankel determinants 
(2.8) A„ := det\pj + k-2]j,k=i,...,n,n > 1, A := 1, 

and 



(2.9) 



£„ := det 



M2n-3 M2n-1 



,n>l, S := 0, 
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defined in terms of the moments {/i ra }n=o,i,...,oo of the weight, 

(2.10) n n := J dx w{x)x n . 
We have integral representations for A n 

(2.11) A„ = — [ dx x ... f dx n f[w(xi) Yl (xk-xjf, n>l, 

' Jl Jl 1 = 1 l<j<k<n 

and £„ 



x k - Xj) 2 , n>l. 



(2.12) £„ = — I d Xl . . . ( dx n f[ w(xi) lf2 x i) II ( 

The three-term recurrence coefficients are related to these determinants by standard 
result in orthogonal polynomial theory |31j , |14j 



(2.13) 
(2.14) 
(2.15) 

with initial values 
(2.16) 



2 A„+iA n _i 
a„ = — , n > 1 , 



->n+l 



A„+i A n ' 
A„ 



n > 0, 



A„+i 



n > 0, 



M0M2 - Mj 
Mo 



7 Mi 2-1 

00 = — , Mo7o = !• 
Mo 



The orthogonal polynomials themselves also have a determinantal representation 

Mo ■ • • M« 



(2.17) 



yj 'A„A n+ ip n (x) = det 



Mn-l 
1 



M2n-1 
2™ 



n > 1, 



and the integral representation 

(2.18) v /A„A„+ip„(a;) = — j dxi . . . j dx n Y[_w{xi)(x—x{) (x k -Xj) 2 . 

n ' Jl Jl 1=1 l<j<k<n 

Another set of polynomial solutions to the three term recurrence relation are 

the associated polynomials {pn\x)}^ =0 , defined by 



n > 0. 



s — X 



(2.19) p^. 1 (x):=jd8w{8)i 
In particular these polynomials satisfy 

(2.20) a n+lP W(x) = (x- b n )p£Ux) - antiUx), 



with the initial conditions p^^x) = 0, Pq~'(x) = Mo7i- Note the shift by one 
decrement in comparison to the three-term recurrence (|2.4j) for the polynomials 
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{Pn(x)}^ =0 - We also need the definition of the moment generating function or 
Stieltjes function 

w(s) 



(2.21) 
(2.22) 



ds 



x£I, 



x — s 

/' 



E^' ^ 7 > x ^°°- 



n=0 



We define non-polynomial associated junctions {e n (^)}$?Lo by 

(2.23) e n (x) :=f(x)p n (x)-p^_ 1 (x), 

which also satisfy the three term recurrence relation (|2.4p . namely 

(2.24) a n+ ie„+i(a;) = (x- b n )e n (x) - a n e n -i{x), 

subject to the initial values e-i(x) = 0, eo(a^) = Jof(x). The associated functions 
have an integral representation analogous to p. 181) 



(2.25) v / A n A„ +1 e„(x) 
1 



(n + 1)! 



dxi 



dx 



w(xi) 



n+1 

n+1 

- LX a; — x; 

;=1 K]<Kn+l 



(Xk-Xj) 2 , X I. 



The polynomials and their associated functions satisfy the Casoratian relation 



p n (x)e n _ 1 (x) - p n -x(x)e n {x) = — , n> 1. 

On 



(2.26) 

Extending (f2~2"]l and (|2~5)l we have 



(2.27) p n (x)=7n 



E b * xn ~ l 



E ^-E«* 

lO<2<J<71 i — 1 



valid for n > 1, while for the associated functions 



(2.28) e„(x)= 7 ,7 1 



Me* 



,i=0 



n+1 



E ^ + E a ? 

iO<i<j<n 



-n— 3 



+ 0(a;- n - 4 ) 



i=l 



valid for n > 0. 



Proposition 2.1 (0,0], [H]). Let 

1 d 



(2.29) 



u>(x) dx 



w(x) 



2V(x) 
W{x) ' 
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for V, W irreducible. The orthogonal polynomials and associated functions satisfy 
a system of coupled first order linear differential equations with respect to x ( ' = 
d/dx) 

(2.30) Wp' n = {fl n -V)p n -a n e nPn - 1 , n>l, 

(2.31) Wy„_i = a„e„_ip„ - (n n + V)pn-i, n > 0, 

for certain coefficient functions V(x),W(x),O n (x),tt n (x). The associated func- 
tions e„,e„_! satisfy precisely the same set of equations. 

If wc define the 2x2 matrix variable 



(2.32) Y n (x;t) 



( v ( x ) e n(x) 

Pn[X) w(x) 

v ,(V) e "-l( x ) 



then the above coupled system can be written as 

/, w d _ 1 / O n (x) - V(x) -a n e n (x) \ vf s 

[ ' dx Yn(x) ~ w(x) { a n e n ^(x) -n n ( x ) - v(x) ) Yn[x) 

It follows that the coefficient functions are specified by 
(2.34) 

®n = W[e n p' n -e' r j> n ]+2Ve n p n , n > 0, 0_i = 0, 
(2.35) 

fi n = a n W [e„-ip'„ - e' n p n -i] + a„V [e n p n -\ + e„_tP„] , n > 1, O = 0. 

Proposition 2.2 (|24j). The coefficient functions arising in Provo sition l 2. 1\ satisfy 
the recurrence relations 

(2.36) {n„ + i-n n ){x-b n ) = W + a 2 n+i e n+1 -a 2 n e n _ 1 , n>0, 

(2.37) 0„+i + Q„ = (x - b n )Q n , n>0 

We will find it necessary to study the zeros of the orthogonal polynomial p n (x) 
which we denote {x\ n < . . . < Xj >n < . . . < x n _ n }. They have an electrostatic 
interpretation as the equilibrium positions of the mobile unit charges, and there is 
a set of equations governing these equilibrium positions known as the Bethe Ansatz 
equations. 

Proposition 2.3 (|16|). The zeros {xj. n }" =1 of the polynomial p n (x) satisfy the 
coupled functional equations 

(2 38) 2 V 1 - 6 »^») W'( Xj , n ) + 2V(x jin ) 

for all 1 < j < n. 
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One can also represent many useful quantities in terms of sums over the ze- 
ros and we illustrate this with an example. Firstly the consecutive ratios of the 
orthogonal polynomials have a partial fraction decomposition 

( 2 39 ) Pn-l(x) = 1 W(Xj, n ) 1 

along with 

(2-40) °» = -£|%T 

Of particular relevance to our application are the semi- classical class of orthog- 
onal polynomial systems [35] defined by the property that V(x) and W(x) in (|2.29|) 
are polynomials in x. The zeros of W(x) define finite singularities of the system 
of ordinary differential equations Q2.30p . (|2.31[) and will feature prominently in this 
study. Let x r be such a point with r > 1. Then at x r the relations (|2.36|) and 
()2.37jl can be combined and integrated to yield 

(2.41) fi*(z r )-V a (s r )=a*e„(x P )e n _i(a: P ) ) n > 1. 

In fact, at a given finite singular point x r , we can deduce the following bi- linear 
identities, that factorise the one above. 

Corollary 2.1. The coefficient functions evaluated at a finite singular point x r are 
related to evaluations of the orthogonal polynomials and associated functions by the 
relations 

(2.42) Q n (x r ) + V{x r ) = 2a n V{x r )p n {x r )e n -x{x r ), n > 1, 

(2.43) Q n (x r ) - V(x r ) = 2a n V{x r )p n ^i(x r )e n (x r ), n > 1, 

(2.44) e„(a- r ) = 2V(x r )p n (x r )e n (x r ). n > 0, 

From the theory of Uvarov the following general result for (|1.29|) is known. 

Proposition 2.4 ( [35 1 ) . The quantity D n (x, x)[w(X)], defined by the equal argu- 
ment form of tl.29\) . is evaluated in terms of the polynomials p n (x) orthogonal with 
respect to w(x) and coefficients j n , A„ of this system as 

(2.45) D n (x,x)[w(X)} = -^—\p n (x)p' n+1 (x) - p n+1 (x)p' n (x)]. 

Proof. This is a specialisation of Uvarov's general result to the case k = and 
I = 2 where the integral D n [xi, X2)[w(X)] is a Hankel determinant with respect to 
the weight wq,2{x), defined by 

(2.46) w ^{x)dx — dp ^{x), Po,2{x) — / (s — Xi)(s — X2)w(s)ds. 
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The non-confluent form of the corresponding identity states that 

(2.47) D n (xi,x 2 ) [w{\) = (A n+2 A n ) ' '-, 

x 2 - Xl 

and the result follows under the confluence x<x — ► x\. □ 

We see from (|2.45|) that our main task is to obtain appropriate characterisations 
of the orthogonal polynomials and their derivatives associated with the weight 

(EH. 

2.2. Deformed Laguerre Orthogonal Polynomials. As we noted in the intro- 
duction we see the appearance of a deformed Laguerre weight (|1.32|) which is a 
member of the semi-classical class with the polynomials V, W in (|2.29|) specified by 

(2.48) 2V(x; t) = -x 2 + (a + 2 - t)x + 2t, W(x; t) = x(x + 1), 

and has finite singularities at x = 0,—t. The moments have the simple evaluation 

(2.49) fi n (t) =t a+n+3 T(n + 3)U(n + 3,a + n + A;t), n > 0, |arg(t)| < tt, 

where U(a, 7; z) is the confluent hypergeometric that is not analytic at z = 0. The 
moments can be written as a sum of two parts one of which is analytic and the 
other non-analytic about t = 

(2.50) finit) = T(a + n + 3) 1 F 1 (-a;-a-n-2;t) 

+ (-ir 3 ^j^ t°^ 3 1 f 1 (n + 3; a + n + 4; t), 
T(a + n + 4) 

where we have to exclude the cases a € 1>q. 

Within the semi-classical class the coefficient functions n (x), f2„(x) are poly- 
nomials with degree fixed independently of the index n. In particular we can relate 
these polynomials to the coefficients of the orthogonal polynomials themselves. 

Proposition 2.5. The coefficient junctions are 

(2.51) ®n(x) = 2n + a + 3 - t - b n - x, n > 0, 

(2.52) Q n (x) = ~lx 2 + l{2n + a + 2-t)x + (n + l)t-a 2 n - 2^1 , n > 1. 

2 2 7„ 

Proof. From the theory of [23] we note that the degrees of the coefficient functions 
are deg0 n < maxjdcgVF — 2, AcgV — 1} = 1 and degi!„ < maxjdcgW — 1, degV — 
1, deg6„ — 1, degt/+l} = 2. To obtain explicit forms for these we use the definitions 
(|2.34jl and (|2.35p and the large x — > 00 expansions of the polynomials and associated 
functions given in ([2~27]) and ([2~2"8]) . The first equalities in (|2~5T|) and ([232]) then 
follow. □ 
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We will also find it convenient to make the following definitions motivated by 
the above result, 



(2.53) 
(2.54) 



2n + a + 3 — t — b n , 
(n + l)t-a n - —. 

In 



Proposition 2.6. The spectral derivatives of the polynomials p n (x) are ( ' 
d/dx) 

(2.55) x(x + t)p' n = (nx + n n - t)p n - a n (B n - x)p n -i 

(2.56) x(x + = a n (0„-i - x)p n - (-x 2 + (n + a + 2 - t)x + t + K„)p„_i 



Proof. This follows from the general form of the spectral derivatives (|2.30[) and 
(|2.30[) . along with the explicit particular forms (|2.51|) and (|2.52p . □ 

Proposition 2.7. The deformation derivatives of the orthogonal polynomials are 
( ■ = d/dt) 



(2.57) 

t(x + t)p n 

(2.58) 

t(x + t)p n -l = —a n (O n -l + t)p n + 



(n + l)t-K n - -(x + t){6 n + t) 



Pn + a n {6 n + t)Pr, 



-(x + t)(e n -! +t) + K n -(n + a+ l)t 



Pn-X 



Proof. We will opt to establish this relation directly from the orthonormality con- 
ditions on the polynomials 



dx w(x)p n (x)p n -i(x) = 8i, , < i < n. 



Differentiating this with respect to t leaves us with the relation 

(2.59) = a / dx w(x) ^ n ^ n ' + dx w(x)p n p n ~i + dx w(x)p n p n ~i, 
J i x + t J j J j 

where use of the logarithmic derivative of w(x) has been made. Now we employ 



Pn- 



-Pn-i + n„_j_i, 



to write the last term of (|2.59p as 



/ dx w(x)p n pn^ l = ^— -5 it0 = I dx w(x)p n p n _i. 

JI Hn—i 



In Jl 
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Considering the first term of (|2.59|) we note 

p n (x)p n -i(x) 



J dx w(a 



x + t 



dx w(x)p n (x) 



Pn-i(x) -p„-i(-t) 



+ Pn-i(-t) / dx w(x) 



x + t 

Pn (x) 



X + t 



-p n -i(-t) J dx w{x) 
= -p n -i(-t)e n (-t). 



^ Pnjx) 
X+t' 



But we can recast p n -i(—t) as 



Pn-i(-t) =^2pn-j(-t)6i,j, 



3=0 



-t) J dx w{x)p n -j(x)p n -i{x), 

dx w(x)p n -j(x) y^^p n -j(-t)p n -j(x), 
3=0 

dx w(x)p n -i(x)^2pj(-t)pj(x). 
3=0 

Combining these we deduce that 

/ dx w{x)p n _i(x) < p n (x) + —p n (x) - ae n (-t)^2pj(-t)pj(x) > = 0, 

for i = 0, . . . , n. The factor in curly brackets in the integrand must be a polynomial 
in x with degree less than or equal to n, and yet is orthogonal to all polynomials 
Pj for j = 0, . . . , n, and thus must be identically zero. This gives us our first form 
for the deformation derivative of p n , 



(2.60) 



Pn{x) = ~—p n {x) + ae n (-t)y^ j p j (-t)p j (x). 

ln 3=0 



Equating coefficients of p n (x) in this relation we find 



(2.61) 



2— = ap n (-t)e n (-t). 

In 



Furthermore using the Christoffel-Darboux formula (|2.7p . and the above equation, 
we can express this derivative solely as a linear combination of p n and p n -i 



(2.62) p n (x) = ae n (-t) 



1 / .s , Pn-i(-t) 



Pn(x) 



I .-sPn(-t) , , 

aa n e n (~t) — p„_i(x). 

x + t 
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Using the bilinear product relations (|2.43[1 and (|2.44|) we arrive at the result (|2.57l) . 
The second of the two relations can be found by shifting in (|2.62| and 

using the three term recurrence relation. □ 

In the matrix formulation the spectral derivative take the particular form 
(2.63) d x Y n (x;t) = (a>o + — + -^-} Y n (x; t) . 



x x + t 

The residue matrices are explicitly given by 

(2.64) .4, = - an9 A, xo = 0,-2 

t \a n V n -l —K n -tj 

(2.65) A t = ± ( {H + 1)4 " *» °f n + t} N V X t = 0, -a 

t \ -a n (6 n -i +t) K n - (n + a+l)t I 



(2.66) 4 



'0 

.0 1, 



Our linear system has two regular singularities at x — 0, —t and an irregular sin- 
gularity at x = oo with Poincare index 1. In this formulation the deformation 
derivative is 

(2.67) d t Y n (x;t) = {l3+^-}Y n (x;t) 

x + 1 I 



(2.68) B=i [-*»-' ° 

2* ^ B n -! + tJ 

As we will see in Section 13.21 p.63| and (|2.67|) form the monodromy preserving 
system corresponding to the fifth Painleve equation with a form equivalent to that 
discussed by Jimbo [H] , in contrast to other forms studied in [T^] , [7] , [5] or [T7] . 

Corollary 2.2. The polynomial coefficients satisfy the following coupled, first order 
mixed deformation derivative and difference equations 

(2.69) 2t— = 2 + b n - X - b n , n > 1, 

a n 

(2.70) tb n =a 2 n -a 2 n+1 +b n , n>0, 

with the initial t = values for b n and a\ given by and \2.88\) respectively. 

This system of differential equations is equivalent to the Schlesinger equations. 

Proof. There are several methods of proof available here. The first is using the 
general result of [24j . expressing the deformation derivatives of the polynomial 
coefficients in terms of a sum of the coefficient functions over the movable finite 



16 



PETER J. FORRESTER AND NICHOLAS S. WITTE 



singular points 



(2.T1) ^ = i£ 



6„(x r ) - e n _i(a; r ) . 



W'(x r ) 



(2.72) 6„ = ^ i?7 ^ x r n>l. 

The only finite singular point contributing here is x — —t. 

Alternatively one can find these derivatives from the polynomial derivatives by 
examining selected coefficients. For example by considering the coefficients of a;™ -1 
in (|2.60[) we have 

(2.73) 7„,1 = ^aen(-*)Pn(-*)7n,i + a^n(-t)Pn-i(-t)7n-i, 
and therefore 

(2.74) (^) = ae n (-iK_!(-i). 

\ In J 

Consequently, using (|2.5p . we find that 

(2.75) — = -a [e„_i(-t)p„_!(-t) - e„(-t)p„(-i)] , 
a„ 2 

(2.76) 6„ = a [a n e„(-t)p„_i(-t) - a„+ie„+i(-t)p„(-i)] , 

which are identical to (|2.69[) and (|2. 70[) respectively. □ 
Corollary 2.3. The recurrences for the polynomial coefficients are 

(2.77) a 2 n+2 -a 2 n 

= 2t + b n+ i [27i + a + 6 - t - b n+ i] - b n [2n + a + 2 - t - b n ] , n > 1. 

and 

(2.78) a 2 n+l [2n + a + 5-t-b n - b n+1 ] - a 2 n [2n + a + 1 - t - 6„_ x - b n ] 

= -b n (b n + t), n > 1. 

The initial data bo(t) and a 2 (t) are given by h2.16\) with the evaluation of the mo- 
ments j2.49\ ). 

Proof. The first relation (|2.77p follows by substituting the explicit forms for the co- 
efficient functions, (|2.5f [) and (|2.52j) . into the relation (|2.37p and requiring equality 
of the polynomials in x. Equality is trivial for x 2 and x , whilst the non-trivial 
equality for a; gives (|2.77j) . The second relation (|2.78p follows from the same pro- 
cedure applied to the recurrence relation (|2.36|) and again the only nontrivial result 
occurs for the x° part. □ 
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This result can also be recovered from a specialisation of work in |5| . In their sys- 
tem of monic orthogonal polynomials the three term recurrence coefficients (3 n ,"fn 
(not to be confused with our use of these symbols subsequently) are related to ours 

by 

(2.79) Pn = b n , j n = a 2 n . 
Equation (39) of this work implies 

n-l 

(2.80) 7n+ i + 7n = (2n + 3)i + (2n + a + 4 - t)(3 n + 2 /?* - 

k=0 

and differencing this once leads directly to (|2.77[) . In addition their equation (40) 
implies 

n n 

(2.81) 7n+i/3n+i = (2n + a + 5 - t - (3 n ) ln+1 + 2Y,lk + MPk + t). 

k=l k=0 

Again differencing this once one finds precisely (|2.78D . 

A check of the above results can be made when t — > whilst all parameters 
are kept fixed since this, as noted before, corresponds to the Laguerre weight with 
parameter a + 2. Therefore we have 

(2.82) al(0) = n(n + a + 2), 

(2.83) b n (0) = 2n + a + 3, 

j!r(j + a + 2) 

(2.84) A n (0)= 3 , 

n! 

In our normalisation we have 

(2.85) p n (x;0) = (-l) w ( = 7 ^ — 1 1/2 4° +2) (*), 

^ r(n + a + 3) J 

where Ln" (x) are the standard associated Laguerre polynomials of degree n and 
index a. We see that the spectral derivative equations (|2.55|) and (|2.56[) reduce 
to the standard expressions for the derivative of the Laguerre polynomials, the 
coefficients in the right-hand sides of the deformation derivatives (|2.57[) and (|2.58|) 
vanish, the recurrence relations (|2.77p and (|2.78[) are identically satisfied, and the 
right-hand sides of (|2.69[) and (|2 . 70[) are zero. 

In fact we will need to develop expansions about t = in order to charac- 
terise our quantities as particular solutions of difference and differential equations 
in Section 3. To this end we have the following result. 
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Proposition 2.8. For fixed n, a ^ Z>o the Hankel determinant K2. <$)) with the 
weight hi. 3 2}) has the expansion about t = 



(2.86) A„ft)=A n (0){ 

1 a /(a- + 1) (a + l)(n-l) 



a 

1 H -ni 

a + 2 



, nt 2 + Oft 3 ) 

4a + 2 V a+1 a + 3 

2I > +1 ) r(a + n + 3) +3 ft 2a+6 )l 

r(a + 3)r2(a + 4) r(n) 1 + Wj + 1 j J' 

mi/i |argft)| < ir. Consequently, under the same conditions, the three-term recur- 
rence coefficients have the expansions about t = 

, , a 2a(2n + a + 3) , „, 

2.87 b n (t =2n + a + 3-— -t+ - ; t 2 + Oft 3 

a + 2 (a + 3)(a + 2) (a + 1) 



r(a + n + 3) fn+3 (1 + Q(t)) + 0(t 2a +6)j 



(a + 2)(a + l)r 2 (a + 3) T(n+1) 
and 

(2.88) a 2 ft) = n(n + a + 2){l - ^ t 2 + Oft 3 ) 

I (a + 3)(a + 2)^(a + 1) 

_ 2 r(a + n + 2) +3 (1 + Q(t)) + ft 2 -+6)l. 

(a + l)r(a + 3)r(a + 4) T(n+1) V y " 1 ; J 

Proof. We adopt the method of expanding the Hankel determinant by expanding 
the moments to leading order 

(2.89) /Lt n (t) = r(o + n + 3) + aT(a + n + 2)t + ^a(a - l)T(a + n + l)t 2 + Oft 3 ) 

, / ^»+i r ( a + l)r(» + 3) 4 „ +a+3 ,„ +a +4\ 

+ l j r(a + n + 4) +Ul j ' 

as t — > using (|2.50p . The determinant can be expanded to leading orders in t, t a+3 
and the resulting determinants evaluated using the identity |26j 



(2.90) det(r(z fc +.?)),- fe=0 ,.. .,„_! = []r(z,) JJ (z fe - 

j=0 0<j<fc<n-l 



n-1 

— TT 



where {zo, . . . , z„_i} is an arbitrary sequence not necessarily in arithmetic progres- 
sion. □ 



We conclude this section by noting some identities relating the polynomial coeffi- 
cients and the zeros of the polynomials. Firstly we give the Bcthe Ansatz equations 
for the zeros of the deformed Lagucrrc orthogonal polynomials which can be directly 
deduced from Proposition 12.31 
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Corollary 2.4. The zeros Xj. n of the deformed Laguerre orthogonal polynomials 
p n (x) satisfy the functional equations 

3 a+ 1 1 „^r- 1 



(2.91) 



u J,n 



t 



3> n 



Or, 



= 1, 1 < J < n. 



According to the electrostatic interpretation, the terms of (|2.91|) can be inter- 
preted in the following way - the first is the interaction of the mobile unit charge 
at Xj tn with the fixed charge of size 3 at the singularity x — 0, the second with the 
fixed charge of size a + 1 at the singularity x = —t, the third with a fixed charge 
of size —1 at the apparent singularity x — 9 n , the fourth the mutual repulsion 
with the other mobile charges and the term on the right-hand side is the linear 
confining potential. From the partial fraction decomposition (|2.39|) specialised to 
the arguments x = 0, — t we have the summation identities. 

Proposition 2.9. The following summations over the zeros have the explicit eval- 
uations 

K n — (n + l)t 



(2.92) 
(2.93) 
(2.94) 
(2.95) 



9„+t 



3i n 

%i ^3,71 



9 n +t 



E 

n 

E 

n 

E 



In addition we can characterise the motion of the zeros with respect to the 
deformation variable. 

Proposition 2.10. The zeros Xj_ n {t) satisfy the differential equation with respect 
to t 

n + t 



(2.96) 



tXj n 



Proof. This follows by equating 

Pn(x) 



(2.97) 



.(*) 



In _ 
In X - X 



and (|23T)) . and then employing (|2~3l))) along with l[2"~3g]) . ([231]) . 

3. Difference and Differential Equations 



□ 



3.1. Difference Equations. In the first subsection we derive an alternative dif- 
ference system in terms of the new variables 8 n (t),K n (t) as specified by (|2.53[) . 
(12-51) . 
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Proposition 3.1. The auxiliary junctions n (t), n n {t) satisfy a system of coupled 
first order recurrence relations 

(3.1) Kn+i + Kn = n (6 n + t — 2n — a — 3), n > 0, 

(3 2) &n 6n ~ 1 = fan ~ t)(K n + t) 

6" Tl 1 6*^ ! + i [k„- (n + a + l)t][K n - (n + l)i]' "~ 
The initial values 9g and Kq are given by 

m) 0(f)- . Jo^dxe^xjt + xy 

(3 ' 3) ^-^f-dxe-xxHt + xy 

Proof. The first of the recurrence relations (|2.77p can be exactly summed and the 
result is 

(3.4) a 2 n+l + a 2 n = (2n + 3)t + (2n + a + 4 - t)b n + 2^6*- b\. 

i=0 

Recalling the second relation of (|2.5p the summation appearing here can done by 
recasting the equation in terms of the new variables and yields 

(3.5) n n+ i + K n = -9 n b„, 

which is (|3.ip . The second member of the coupled set is most easily found from 
the general relation (|2.41[) evaluated at the finite singular points x = 0, — t and 
employing the new variables. These two key identities are 

(3.6) (K n + t)(n n - t) = a n 6 n 6 n -i, 

(3.7) [ Kn -{n + a+ l)t\[n n - (n + 1)<] = a 2 n (6 n + + t). 

The ratio of these two identities yields the relation (|3.2[) . □ 

There are other recurrence relations which will be used subsequently, and the 
first is 

(3.8) a n {6„ + 0„_i +t) = -(2n + a + 2)k„ + [n 2 + (n + l)(a + 2)}t. 
This follows from the subtraction of (|3.6p from (|3 . 7[) . The second relation 

(3.9) a 2 n+1 - a 2 n -b n -t = 2n n + b n 6 ni 

is derived by writing the definition of — n n in terms of the old variables and 
then employing (|3.ip . The last relation 

(3.10) a 2 l+l 6 n+ i - a 2 n 6 n -i = b n (2n n + b n 9 n ), 
is a consequence of (|2.78p along with the use of (|3.9[) . 
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As a consequence of relations (|3.6[) from (13. 7p we have 



Kn. t 



9 n + t 



n ™n 



t K n — (n + l)t 



2 "n-l 



a 



On + t 
2 ®n-\ + t 



t \ n K n + t n n n - (n + a + l)t J ' 



(3.11) 

and 

(3.12) 



t 



K n t 



K n — (n + a + l)t 

K n + t 



n + a + 2 + 



K r , + t 



+ { n + a + 2 + 



>n-\ 



9 n + 1. 



m— >n+l 



3.2. Reduction to Painleve V. Here we will identify the fifth Painleve system 
as the solution to our system of equations characterising the deformed Laguerre or- 
thogonal polynomial system. This is most simply seen in terms of the new variables 
6 n , Kn rather than the basic orthogonal polynomial variables a n ,b n . 

Proposition 3.2. The auxiliary quantities n (t), K n (t) satisfy the coupled first or- 
der ordinary differential equations 



(3.13) 
and 

(3.14) tk Tl 



? n +t 



tO n = 2n n + 9 n (2n + a + 3-t- On), 



+ -j- }K 2 n + \2n + a + 3 - (2n + a + 2)- 



- [n 2 + (n + l)(a + 2)]t - — + (n + l)(n + a + 1) 



+ t 



Equations i3.13\) and {3. 1J$ can be solved in terms of the fifth Painleve system 

<1 



(3.15) 



t- 



K n = t(l + qp), 



where q,p are the Hamiltonian variables of the Okamoto Py |27j system with the 
parameters 

„2 1 



a = y, p = -2, 7 = -(2n + a + 3), 6- 



(3.16) 
or 

(3.17) V2 — v\ = —2, V3 — v\ = n + a, V4 — v\ = n, V3 — v± = a. 

The solutions satisfy the boundary value data at t = 
2 



(3.18) n {t) 



2a(2n + a + 3) ^ + 



t->0 a + 2 (a + 3)(a + 2) 2 (a + 1) 



(a + 2)(a + l)r 2 (a + 3) r(n+l) 



r(a + n + 3) ta+3 (1 + Q(t)) + 0(t 2a+ 6)) 
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(III 



d 



(3.19) Ml=^t + + \+?° * + 0<!*) 

t-*o a + 2 (a + 6) [a + 2) z ya + 1) 



r(0 + n + 3) fa+3 (1 + Q + 0(i 2a +6)i 



(a + 2)(a + l)r 2 (a + 3) T(n) 
provided a (£ Z>o and |arg(t)| < tt. 

Proof. If we employ (|3 . 9() in (|2.70p and then substitute for b„ in terms of 9 n then 
the result is (|3.13|) . Let us define the shorthand notation 

(3.20) r„ := 

In 

Furthermore when the deformation derivative (|2.70|) is summed on the free index 
the result is 

(3.21) tt n =al + T n . 

Now if compute the deformation derivative of K n and use (|2.69j) along with (|3.21|) 
we arrive at 

(3.22) tk n = K n - a„{9 n - 0„_i). 

Now the idea is to eliminate a n and 0„_i, which appear in (|3.22[) . through use 
of the recurrence relations. Equation (|3 - 2|) is a linear equation for & n -\ in terms 
of the unshifted variables and the solution can be substituted into ()3.8j) yielding 
a linear relation for a n in terms of unshifted variables. Then both solutions can 
be substituted into (|3.22j) and the result is ()3.14|) . One can easily verify that the 
transformation to the Hamiltonian variables q,p (|3.15p with the parameters (|3.16|) 
yields the Hamilton equations of motion for the Hamiltonian in |27j . □ 

Remark 3.1. A few remarks can now be made regarding the identification of the 
recurrences (|3.1[) and (|3 . 2|) . This is different in appearance from the discrete in- 
tcgrable equations that arose in the study of the Laguerre unitary ensemble |13j 
which were explicitly identified with the system in the Sakai scheme, with the ra- 
tional surface — ► and has a continuous limit of Painleve IV . In fact we 
find that the variables of the latter system can be expressed in terms of our own 

/0 r, \ K„ + 9 n (n + a+l-t-8 n ) t 

(3-23) x n = — — , y n = -—, 

t>„ + 1 u„ 

and it is clear that one cannot transform (|3.ip and (|3.2[) into this system using 
such a transformation. Also the two systems arise as different Schlcsingcr-type 
transformations - in our case as a sequence where ao t— > ao + 1,0:2 ^ 0:2 — 1 
whereas in the other case as oq >-* 00 + 1 , 03 i— > 03 — 1 . 
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Remark 3.2. In [13j two fundamental quantities were studied - the r-function r[n](t) 
and its logarithmic derivative the cr-function V n (t;a, fx). Their relation to the ob- 
jects of the present work are 

(3.24) r[n] = c(n, a)n! e -™*t- n2 -™ (a+4) A n (t), 
where c(n, a) is an unspecified constant and 

(3.25) V n (t;a,2) = -nt - An + t^ log A n (t). 
We also note that 

(3.26) V n (t; a, 2) = T n + n(n + a-2-t) 
and consequently 

(3.27) r„ = -n(n + a + 2) +t^- log A n (t) 

at 

(3.28) b n = 2n + a + 3 + t^-\og 

at /\ n+ i(t) 

Remark 3.3. The new variable T n (t) possess an expansion as t — > which can be 
directly found from ([2~86]) and ([237]) 

„„v / ^ , % a 2n(n + a + 2) a 9 q. 

3.29 n(n + a + 2 + T n (t) = —-nt - ——±——±—t 2 + O t 3 

a + 2 (a + 3){a + 2)- (a + 1) 



r(a + n + 3) ^ a+3 



(a + 2)(a + l)r(a + 3)r(a + 4) T(n) 
again provided a (£ Z> . 



(l + 0(<)) + 0(i iQ+b ), 



Remark 3.4. The use of Proposition 13.21 is in the computation of the orthogonal 
polynomials in (|2.4|) corresponding to the weight (|1.32p . For this we note from 
(|2~53)) that 

(3.30) b n = 2n + a + 3 - t - 6 n , 
while (|23i|) together with (|23| . (|23|) show 

(3.31) a? = (n + l)t- — - K n 

In 
n-1 

(3.32) = (n + l)t + Y^h ~ K n , 

j=o 

(the quantity a n is positive, so the positive square root of this equation is to be 
taken). Further, it follows from (|2. 13[) and (|2.5p that 

(3.33) ^ o o = A nj 

7o7r---7n-i 

where each 7j is the coefficient of x J in as specified by (|2.2j) . All terms in the 

equation (|2.45|> for D n {x,x) are then known, and the task is then to compute the 
integral as required by (11.311) . 
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An alternative system of coupled first order ordinary differential equations which 
will be used for scaling to the hard edge is given in the following proposition. 

Proposition 3.3. The variables 9 n (t),T n (t) satisfy the coupled first order ordinary 
differential equations 

(3.34) W n = 6 n - \ei - 2(2n + a + 2 - t)6% 

+ [4T„ + (2n + a + 2 - t) 2 - 4(n + l)t]6 2 n 

+ 4t[r„ + n(n + a + 2-t)+a + 2- t]d n + 4t 2 | 



1/2 



(3.35) tt n = (n + l)t + -e n (2n + a + 2 - t - 6 n ) + - 2(2n + a + 2 - 
+ [4r„ + (2n + a + 2 - i) 2 - 4(n + l)t]6i* 

+ 4t[r„ + n(n + a + 2-t)+a + 2- i]6>„ + 4t 2 [ . 

Proof. We proceed in a series of steps. Firstly we use (|3.6[) to solve for in 
terms of n , k„ and T„. In the second step we substitute this solution for 9 n _i into 
(|3.7p and solve the following quadratic equation for n n in terms of 9 n and r„, 
(3.36) 

Kl+e n (2n+a+2-t-9 n )K n + [(n+l)t-T n }9 n (9 n +t)-[n 2 + {n+l)(a+2)}t9 n -t 2 = 0. 

The choice of the sign of the square-root branch follows from the expansions (|3.18|) 
and (|3.29p on one hand, and on the other hand noting that as t — > 

(3.37) \9i-2{2n + a + 2-t)9l 

+ [4r„ + {2n + a + 2-t) 2 - 4(n + l)t]6> 2 

1/2 



4i[r„ + n(n + a + 2-i) + a + 2- t]0 n + 4* 



2 



2a(2n + a + 3) , ^ 2 + Q(t 3 ). 



(a + 3)(a + 2) 2 (a + l) 
In the final step we use these solutions for 9 n -\ and K n in (|3.22|) and (|3.2ip . □ 

In preparation for the hard edge scaling limit we need to make evaluations of 
the polynomials at the finite singular points. Firstly considering x = we note 
that 

,ooo> n n(0) PrJfi) Kn + t 2 n 

(3.38 J — - = a n — = — = a n -, 

7r„_i(0j Pn-liP) Qn-1 Kn — t 

as follows immediately from (|2.55|) . Furthermore we also have 
(3.39) t(lo g7 r„(0)) = i(logp„(0)) + he n + t)=n + 

Z tin. 
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which follows from (12.571) and the above result. The corresponding result for the 
polynomial ratio at x = — t is 



(3.40) 



p n (-t) 1 K„-(n + a+l)t_ 6 n + t 



(3.41) 



Pn-l(—t) «n On-l + t K n - (n + l)t 

After [T] we define the orthogonal polynomial ratios 

Pn(x;t) 



)n(x;t) 



Pn(0;t)' 

because we are interested in the scaling properties of the orthogonal polynomial 
system at the edge of their interval of orthogonality, x = 0. Then (|2.4[) and Propo- 
sitions [2J3 and [277] can be translated into the following three corollaries. 

Corollary 3.1. The three-term recurrence for {Q n }n=o,i,... system is 

(3.42) b n (Q n+1 + Q„_i - 2Q n ) + (b n + 2^—±)(Q n+1 - Q n _{) + 2xQ n = 0. 

Corollary 3.2. The spectral derivatives of Q n ,Q n -i are 
(3.43) 

X 

x(x + t)Q' n = nxQ n + («„ - t) Q n - Q n -i + -r-Q n - 
(3.44) 

x(x + t)Q' n _ 1 = x[x - (n + a + 2 - t)\Q n -x + («n + t) Q n — Qn-i 



Corollary 3.3. The deformation derivatives of Q n ,Q n — i are 

(3.45) t(x + t)Q n = -x(n + ^-^)Q n + (k„ - t)^-^ [Q„_i - Q n ] , 
(3.46) 

i(x + t)Q n -x = x[n + a + 2 + K " + - ]Q n -i + (k„ + t) ^"^ 1 + - [Q n -i - Q n ] ■ 
As we noted earlier the polynomial ratio Q n (x;t) has a product representation 
(3.47) 



Qn(x;t) = H[l- 

3=1 



where again Xj <n is the j-th zero of the polynomial. We can use this fact to compute 
sums of the inverse powers of the zeros from the above differential equations. 

Proposition 3.4. The increment of the sum of the reciprocals of the zeros going 
from n — 1 to n is given by 
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Proof. The required increment of the sum of the reciprocals of the zeros is the 
order x term in the expansion of tp n {x; t) := Q n {x] t)/Q n -i(x; t) about x = 0, and 
this can evaluated from the same expansion of the differential equation for ip n with 
respect to x. This latter differential equation is easily found from (|3.43|) and (|3.44[) 
and is 
(3.49) 



x(x+t)ip' n = - (K n -t) + (2K n + x[-x+2n+a+2-t])tp n + 



x - (/„_! 



(k„+*)^. 



□ 



At the other finite singular point, x = —t, we have as a consequence of these 
corollaries 

Q n (-t;t) K n ~t n + t n n -(n + a+l)t 9 n -i 



(3.50) 
and 
(3.51) 



Q n -i(-t;i) K n -(n + l)t 



*-!+«' 



-Q n {-t-t) = -^——JtQ n (-t;t). 



Note that the derivative with respect to t in the latter equation is a total derivative. 

To complete our preparations for the hard edge scaling we need to identify two 
polynomial variables that will scale to independent variables in the scaling limit. 
The first is the orthogonal polynomial ratio Q n , and for the second a number of 
choices could be made but a simple choice is 



(3.52) 



R)i ■ — Qn Qn—X* 



Corollary 3.4. The spectral derivatives ofQ n ,R n are 

K n t 



(3.53) x(x + t)Q' n = x [ n + 

(3.54) x(x + t)R' n =x 



Qn ~t~ \fi>n -^Ti; 



n -\ h n + a + 2 H — x — t 



Q, 













0n ; 



+ x(a; + 7j) - (a + 2)x-2t 



Proof. This follows from Corollary! 

Corollary 3.5. The deformation derivatives ofQ n ,R n are 
(3.55) 

t(x + t)Q n = -X \Tl + ^ -J Qn - («n - t)-^7r—Rn, 



□ 



(3.56) 

+ <)-R,i = —x 



+ ?i + a + 2 + -4 



a; I n + a + 2 



Km + I 



+ (Kn+t)- 



(«n - t)- 



+ t 
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Proof. This follows from Corollary 13. 31 □ 

3.3. Inequalities and Bounds. A key step in proving our hard edge scaling limits 
will be bounds on the variables 6 n , n n and some auxiliary quantities. The first step 
is the following result. 

Lemma 3.1. The variables O n (t),K n (t) satisfy the inequalities 

(3.57) 6n r +t - < -A_ < 0. 

K n — {n + l)t K n — t 

(3.58) AzL < 9 r 1+t - < 0, 

K n + 1 K n — (n + a+ 1)1 

for all positive, real and bounded t and n > 1. 

Proof. That the ratios given in (|3.57p and (|3.58|) are negative is a consequence of 
the fact that the polynomial p n (x) evaluated on the negative real axis, i.e. exterior 
to the interval of orthogonality, has a fixed sign. Specifically (— l) n p n (— y) > 
for real, positive y. Using the ratio relations (|3.38[) and (|3.40[) we have the upper 
bounds. From the Christoffcl-Darboux formula (|2.7p at equal arguments we note 
that 

(3.59) Pn-^p'nix) ~ p' n - 1 {x)p n {x) > 0, 

and from the above p n -i{x)p n {x) < for x G — R+ we conclude 

(3.60) < ^44, 

Pn{X) Pn-l(X) 

under the conditions on x. Integrating this inequality from to — y 6 — R + we 
arrive at 

/o £i \ Pn{-y) ^ Pn(0) ^ „ 

(3.61) — - < -— < 0. 

Pn-l(-y) Pn-l(0) 

Then identifying these ratios with p.38p and (|3.40[) in the case y = t leads to the 
relative inequalities. □ 

The above set of inequalities must all apply simultaneously and we see in fact 
that it implies restrictions on the variables 9 n , n n . 

Lemma 3.2. For bounded t £ R + and n > the variables 9 n , k„ satisfy inequalities 
which place them in one of three cases, as illustrated in Figure 1 - 
Case I: 

(3.62) < e n , 

(3.63) 9 n {9„ +t-2n- a - 3) +t < K n <t- n6 n , 
Case II: 

(3.64) -t < 9 n < 0, 

(3.65) t - n0 n <n n <t + min{nt, 6 n (6„ + 1 - 2n - a - 3)} < (n + l)t, 



2>s 




-t 



Figure 1 . A pictorial form of the inequalities taking the example 
of a = 1/2, n = 2 and t = 5/3, which illustrates the generic 
situation for n + a + 3 > t. 



Case III: 

(3.66) 9 n < -t, 

(3.67) (n+ l)t < k„ < t-n6 n . 

Proof. From the three inequalities implied by (|3.57|) we see that 9 n ^ according 
as K„ ^ t, 9 n + 1 ^ according as K„ ^ (n + l)t, and 9„(9 n + 1) ^ according as 
Kn ^ t — nO n . In (|3.58p we make the replacement m-m+l and employ (|3.ip to 
eliminate This inequality now reads 

(3.68) - t - v ; < ^ < b„ = 2n + a + 'i-t- 9 n . 

Consequently these three inequalities imply 9 n ^ according as n n — t ^ 9 n (9 n + 
t — 2n — a — 3), n + 1 ^ according as n n — (n + \)t ^ (#„ + t)(#„ — 2/i — a — 3), 
and 9 n (9 n + t) ^ according as k„ ^ 9 n (9 n + t — n) + t. Combining these sets of 
inequalities leads to the three cases above. □ 

We see that Case II applies in our situation. 

Lemma 3.3. For all n and t G R + we have —t < 6 n < and t — nO n < K n < 
{n+l)t. 

Proof. From the residue formula (|2.44[) we recall that 0„(O) = 9 n = 2tp n (0)e n (0) 
and O n (— t) = 9 n + t = —atp n (—t)e n (—t). As we noted in the proof of Lemma 
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(|3.ip it is immediate from the integral representations of the polynomials and their 
associated functions, ([2l8]) and ([2725]) . that (-l)> n (-t) > and (-l) n+1 e„(-i) > 
for all real t > 0. This places 9 n in the range applying to Case II. □ 

We can also draw some conclusions concerning the zeros of the orthogonal poly- 
nomials which will be important subsequently. 

Corollary 3.6. Each zero Xj :H (t) is a monotonically decreasing function oft and 
interpolates between the Laguerre zero with t — and exponent a + 2 and the 
Laguerre zero with t = oo and exponent 2, 

(3.69) x i,n(0) > Xj,n(t) > x 3 , n (oo) > 0, 

for all 1 < j <n and bounded t > 0. 

Corollary 3.7. The following bounds on the reciprocal sums over the zeros hold 



x j,n{t) 2 



n 

< 



n 



n 

< 



.r- 



J Xj^ n (t) +t a + 



Proof. From the two-sided bound on n we can deduce 
(3.72) 1 < -\ — — < 



t (t) + t Xj, n {t) - 9 n Xj,n{t) 

Employing this in the Bcthe Ansatz (|2.9ip summed over j we arrive at the above 
bounds. □ 



4. Special Case a e Z> 

Our evaluation of the distribution function in terms of the fifth Painlcvc sys- 
tem is with all three free parameters variable in some sense - one is fixed in this 
application at a positive integer, one is the index n£Z and the remaining one is 
a G C. Up to this point we have studied in some depth the recurrence relations 
with respect to n while a has been left arbitrary other than being restricted because 
of the existence considerations. From the point of view of the Painleve theory it is 
quite natural that the transcendental objects become classical when a G Z for cither 
positive or negative subsets of the integers. In particular it is expected that the r 
functions in the theory will have Hankel dcterminantal forms of classical function 
entries with a rank dependent on a. It is these cases which have been studied in the 
past |12j.|llj using methods which transform the integral into the determinantal 
representations and then employ confluent Vandcrmonde identities. 



:S0 



PETER J. FORRESTER AND NICHOLAS S. WITTE 



Proposition 4.1. When a € Z>o we have the evaluation for the Hankel determi- 
nant 



(4.1) 

and for a = 
(4.2) 



A„(t) 



(n + a)\ 



A„(i) 



Cn+l,n+l 



Proof. In [12] Eq. (3.18) (after correcting) states 

(4.3) A n (t) = ^^gp (-1)^^ det[23J+*- a X^l? (x)!,^]^!,...^ 
where := d/dx. Using the Laguerre polynomial identity 

(4.4) Z3™iW(x) = (-ir4 a _ + r T ) W ! mGZ> 0! 
with the proviso Li? (x) = for n < 0, we arrive at (|4.1j) . 



□ 

As a consequence of the relations (|3.28j) and (|3.27j) the variables 6 n (i), T n (t) will 
have a x a determinant forms, and in particular for a = 

(4.5) e n (t) = -t, n n (t) = (n+l)t, r„(t) = -»(n + 2). 

The orthogonal polynomials also have determinantal forms of the following type. 
Proposition 4.2. When a € Z>o £/*e orthogonal polynomials are given by 

(4.6) v/A„A n+1 p„(x; i) = (-l)«+»+L4 i J a ! ...(„ + a )m ...(„+ i)! 

r (fc+l-a) / 



x (a; + t)" a det 



j-U+k-a) / i \ 



fc=l,...a+l 

j — l : ...,a 
fe=l,...a+l 



(n + 2)(n + l) 



1/2 



and /or a = 

(4.7) p„(x;t) = (-l)" 

If a > Ti + l we note that £„°M_ a (— i) = 0. Consequently, under the same condition, 
the polynomial ratio is given by 



dot 



(4.8) 



i(x;t) 



x + 1 



T (k+l-a) , x 

r (j+k-a) /_.\ 
1J n+a+2-j-k\ L I 



fc=l,...a+l 

J — l,...,a 
fe=l,...o+l 



det 



n+2 
n+a+1— fe 
j(j+k-a) 

^n+a+2-j-ky V 



fe=l,...a+l 



J — l,...,a 
fe=l,...o+l 
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Proof. Starting with the integral representation (|2. 18|) we follow the procedures 
used in [12] . Taking one factor of the squared product of differences we write it like 

( 4 - 9 ) II s * ~ x i) = n n-i de ^4-i(^)kfc=i,..,rn 

\<j<k<n 1 lz=0 Cl 

using the Vandermonde identity and where c n is the leading coefficient of L^"\x) 
and a is a parameter to be fixed later. Of the remaining factors in the integrand 
we write 



(4.10) \{xj+t) a {x-Xj) (x k - Xj ) 



{ ' N -X (» + *)"" 



l<j<k<n 



ULa ll U 1 ilo ± ci 



x det 



-(a) 



j — 1 ,n 

k=l,...,N 



7-1 r( Q ) 



fe=l iV 



j — 1 , a 

fc=l,...,iV 



where the confluent Vandermonde identity has been used and N = n + a + 1. 

Reassembling the integral with these two factors, then expanding the determi- 
nant in (|4.9p we multiply each of n factors into the determinant of (|4. 10[) . Making 
use of the antisymmetry of the row ordering in the first n rows of the determinant 
we can perform the n integrals as long as we choose a = 2. Then 



n-l 



(_l)a(n+l) 

(4.11) y/A n A n+lPn (x; t) = ttq _ 1 - rTjv _ 1 [] 



U + 2)! 



=0 Cl 7=0 



(x + t)- a 



x det 



-(2) 
J n+k 

7-1 r(2 



-l0«0 



fe=l,...,o+l 



^- lL n+ fe -l(2/)l,=-* 



j=l,...,a 
fc=l,...,a+l 



Using the identities 
(4.12) 



L^- 1 \x)=L^{x)-L { :\{x), 



along with elementary column operations, and then identity (|4.4p we are lead to 
(|4.6[) . The evaluation for the polynomial ratio (|4. 8[) is a simple consequence of the 
first evaluation along with 



(4.i3) 4+;+i + -2(o) = 



□ 
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Proposition 4.3. When a £ Z>o the deformed Hankel determinant is given by 
(4.14) D n (x,x) = c n+2 , n+2+a (-l)^ a+1 ^ a+2 Hx + ty 2a 



x dct 



T (J+k-l-a) ,_,x 
^n+a+3-j-k\ u ) 

~ r{j+k-l-a) 
J n+a-\-3—j—k 



(■r) 



j — l....,a 
fe=l,...o+2 

3=1,2 
fe=l,...a+2 



x dct 



Proof. The quantity D n (x,x) was essentially computed in [TJ] in Eq. (3.20), which 
can be recast as 

(4.15) D n (x,x) = c n+2 , n+2+a (-l)^ a+1 K a +V(x + t)- 2a 

Di+ k - 2 L-i a -^(u)\ u=x 

Then (|4.14[) follows by application of the identity (14. 4p . 

As a consequence the distribution of the first eigenvalue spacing is 



(4.16) A n>a {y) = (_l)^+l)(«+2y e W dtt a( t _ y ya e -(n+2)t 

Jy 

T (J+k-a-l) ,_ ,\ 



j — l,...,a 
fc=l,...a+2 

3=1,2 
fc=l,...a+2 



□ 



x dct 



r (j'+fe-a-l) / x 



j — 1 ..... a 
fc=l,...a+2 

3=1,2 
fc=l,...a+2 



for a = 1, 2, 3, . 



and 



(4.17) A nfi (y) 



for a = 0. 



-t/ 2 e _( " +1)y 



5. Hard Edge Scaling 
5.1. General Case. We define new scaled spectral variables s, z by 
(5.1) 



s 

t= 4n"' 



z 

in' 



in the triangular domain s > z > and study the scaling of the finite distribu- 
tion (|1.31|) as the polynomial degree n — > oo. What is required here is not just 
the asymptotic scaling of the orthogonal polynomial coefficients but also of the 
polynomials themselves in the neighbourhood of an endpoint of the interval of or- 
thogonality. For the deformed Laguerre polynomials this would be a generalisation 
of the asymptotics of Hilb's type for the Laguerre polynomials as found in Szego 



(5.2) 



n\{\Myl 
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with M = An + 2/i + 2 as n — ► oo and Iu(z) the standard modified Bessel func- 
tion. Despite a resurgence in activity around these questions, especially the use of 
Riemann-Hilbcrt techniques on these problems, there are no results available for 
our particular problem. A general review of the asymptotics of orthogonal polyno- 
mials can be found in (23) . and an introduction to the Riemann-Hilbcrt approach 
to the asymptotics is the chapter in [22] . However to establish the nature of and 
the existence of limits for our variables we do not require such techniques. 

Lemma 5.1. Under the scaling of the &5.1\ ) An6 n (t) and n n (t) are bounded for all 
real positive t and n > 1 . 



Proof. From the result of Lemma 13.31 we sec that 

(5.3) - s < An6 n (s/An) < 0, s/An - AnO n {s/An) < n n {s/An) < (n + l)s/4ra, 
and the assertion follows. □ 
Corollary 5.1. Under the above conditions 

K n (t)-t 



(5.4) 



= 0(1), 

t=s/4n 



9 n (t) 

Proof. The above lemma states that 

(5.5) 6 n (s/4n) = 0(1), asn^oo, 
and we find as a consequence that also 

(5.6) e := n(n + a + 2) - s/A + T n {s/4n) = 0(1), as n -> oo. 



The discriminant T> appearing in the workings of Proposition ^. 3l can then be written 
as 

(5.7) 

V 2 := 6^-2(2n+a+2-t)6l+ [A(e-t) + {a + 2-t) 2 -Ant}d' 2 l + At[s + a + 2-t}d n +At 2 1 



and therefore 
(5.8) 



V(t) 



= 0(1), as n 



On{t) 

From the formula for K n in Proposition 13.31 we note that 



,r „n K n — t a+2 t 1 , . . V 
(5-9) n + -fT- = —-TT + ^n+t)-—, 

t) n Z tl n Z ZV n 

and the result then follows. □ 

Proposition 5.1. For bounded s G R+ under the scaling 15. 1\) the variables n {t) 
and T n (t) converge to limits in the following manner 

(5.10) lim 4nO n (t)\t=s/4n = Ms), 

n — 'oo 

(5.11) lim n(n + a + 2) + T n (t)\ t=8 / An = v{s). 
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The variable K n (t) converges like 
(5.12) 

Proof. Firstly we note that 



lim «„(t)|t= a /4n = ~tM s )- 



rj := a 2 „(t) - n(n + a + 2)| t=s/4n = 0(1), 



asn^ oo, 



as follows from (|5.6p . Starting with (|3.8p we see that 

(2n + a + 2)n n - [n 2 + (n+ l)(a + 2)]t 



%i — On-1 — 20 n + t 



1 



= <^ (2n + a + 2 + 7?)t + 2(n + a + 2)0„ n + 



as n — > oo. 



(5.13) 
and thus 

n0 n -(n-l)0 n _i = 0„_i+-^- ^ (2n + a + 2 + rj)t + 2{n + a + 2)0„ ( n + 

a n I 

Thus we have shown 

n9 n {t) - (n - l)0n-i(*)U 8/4 „ = 0( ? i" 1 ), 
Now we can write the quantity of interest 

n °n{t)\ t=s/in - (n- l)0n-lW| t=s /4(„-l) 

= [n0„(i)-(n-l)0 n _ 1 (i)]| t=s/4 „ 

+ (n - l)0«-i(*)| t=s/4n - (« - l)^i(i)lf=,/vi) . 

so we require bounds on the difference of the last two terms on the right-hand side. 
Let t = s/An and £> = s/4(n— 1). Because 0„(i) is continuously differentiable and 
its derivative is given by (|3.13D 

(n-l)|0„_ 1 (*)-0„_ 1 (i > )| 

<(n-l)(i>-i) max |0„_i(m)|, 
«e(t,«>) 

< — max l^^^u)!, 
4n ue(t,t>) 

< — max u~ 1 |2k„_i(u) + n _i(it)(2n + a + 1 — u — n _i(it))|, 
4ri ue(t,*>) 



< — max u 1 |0„_i(w)| 
4n ue(t,t>) 



2 re - 1 



«n-i(w) - u 



< i> max 
= 0(n 



«6(t,t>) 
-1 



n - 1 + 



-iW 

«n-l(«) - 



+ a + 3 — u — n _i (u) + 



2« 



>n-l(«) 

+ a + 3 + u+ |0„_i(u)| + 



2u 



|0n-l(u)| 



as n — > oo. 
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Thus the limit shown in (|5.10p exists. Turning our attention to (|5.12|) we can use 
()3.ip to compute that 

^n+l 2iK n b n 9 n , 

= -2n n - (2n + a + 3 - t - 9 n )9 n , 

K n t 



-29„ n + 



2t-(a + 3)9 n + 9 n {9 n + t), 



as n — > oo. 



and therefore 

[n n+ i(t) - «;„(*)] | t=s/4n = 0(n" 1 ), 
Now in this case the quantity we require is 



K "+l(*)lt =s /4(n+l) K n{t)\ t=s / 4n 

= «n+l(*)| t=s / 4 ( n+ l) ~ K «+l( t )lt= s /4n + [ K n+l(t) ~ «n(*)] l t=s /4n ' 

and therefore we need to bound the difference of first two terms on the right- 
hand side. Let us denote t< = s/4(n + 1). Again because K n (t) is continuously 
differentiable with derivative (|3.22|) we have 

\ K n+l (t<) (*)l 

< (i- t<) max 
ue(t<,t) 



< 



4n(n + 1) «e(t<,t) 



max \k n+1 (u)\, 



< 



max u ^/tn+i^) - a,n +1 (u)(9 n+ i(u) - 9 n (u))\, 



An(n + 1) U £(t<,t) 

< — max f [« n+ i(u)| + (2n + a + 4 + |jj(u)|)u 
n «e(t<,t) V 

«n+i(") - u 



+ 2(n + a + 3)|(9„ +1 (V)| 
0(n _1 ), as n — > oo. 



n + 1 



9„+i(u) 



In the last two steps we have used (|5.13[) and the subsequent estimates. Thus the 
limit in (|5 . 1 2|) follows. The fact that the limit of n n (t) under the hard edge scaling 
is related to the limit given in (|5. 10[) follows from the relation (|3.2| . The limit 
(|5.1ip is a consequence of the limits in the primary variables. □ 



In addition the following combinations of variables possess scaling limits which 
will subsequently be useful. 



Mi 
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Corollary 5.2. For bounded s E the following limits as n — » oo exist 

2K n (t) + e n (t)b n (t) 



(5.14) 
(5.15) 
(5.16) 
(5.17) 



n + a + 2 + 



K„ (t) — t 

n I nH 7TT1 \-n + a + 2 + 

a n \t) 



Unjt) - t 
On(t) 

Kn(t)+t 
0n-l(t) 

K n (t) + t 



t=s/4n 



t=s/4n ' 



t=s/4n ' 



t=s/4n ' 



C(s), 

) 

-C(s), 



Proof. The limit in (|5.15[) is a consequence of (|5.4[) in a previous corollary. The 
scaling limit of (|5.16p follows from that of (|5 . 1 5[) and the identity (|3 . 1 1 1) . The limit 
(|5.17p can be derived from (|3 . 1 2[) and as a result one can deduce that 

i Kn ~ t „ K n + t 

n H h n + a + 2 H — 



is of order 0(n 1 ). 

We note other relations amongst the scaling limit functions, namely 

(5.18) 2/j,C(s) = -[(a + 2)fj, + 2s] 

- {[(a + 2)/i + 2s] 2 + ifj,(fj, + s)v - + s) 2 } 1/2 , 

and 

sC(C + a) 



□ 



(5.19) 
and 
(5.20) 



2C + a + 3 = s 



H + s 
ji-2 



I' 



Proposition 5.2. The scaled variables fj,(s),u(s) are characterised by solutions to 
the Pup system with parameters V\ — a + 2, v-2 — a — 2. In particular 



(5.21) 



u{s) = -aiii'(s) + - a - 2, 



where <Tni'(s) satisfies the Jimbo-Miwa-Okamoto a-form for Pnr with above param- 
eters. The boundary conditions to uniquely specify the solution v(s) are 



(5.22) v(s) 



+ 8(a + 3)(a + 2) 2 (a + 1) 

2 

~ 4 a + 3 (a + 2)(a + l)T(a + 3)r(a + 4) 
assuming a (£ Z> and |arg(s)| < 7r. 



s 2 + 0(s 3 ) 

s «+3 (i + o(s)) + 0(s 2a + 6 ), 
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Proof. Introducing the scaling ansatzes (|5.10j) and (|5.1ip formally into the differen- 
tial equation (|3.34p under the scaling (|5.1|) we find that the highest order nontrivial 
relation 

(5.23) /x(s) + s = 4si/(s), 

at order n" 1 . Proceeding in the same manner with the differential equation (|3.35D 
we find the highest order relation is 

(5.24) sfj(s) = n + {[(a + 2)(i + 2s} 2 + 4/i(# + s)v - n(fi + s) 2 } 1/2 , 

which occurs at order n -1 as well. Eliminating fi(s) using (|5.23|) we find that (|5.24|) 
yields 

(5.25) s 2 (i)) 2 - (a + 2) 2 (i>) 2 + i>(4i> - l)(sz> - v) + -a{a + 2)v - —a 2 = 0, 

2 16 

which is almost the Jimbo-Miwa-Okamoto cr-form for Pnr [28]. The boundary 
conditions follow from the application of the scaling limit ([5. lip to the expansion 
about t = 0, Equation ([3~29]) . □ 

In addition we find the scaling behaviour of the Hankel determinants and poly- 
nomial evaluations to be given by the following propositions. 

Proposition 5.3. As n — > oo under the hard edge scaling the Hankel determinants 
scale as 

(5.26) A„(i)| t=s/4 „ ~ l!...(n-l)ir(o + 3)...r(n + a + 2)exp( / —v(u)\, 
and the monic polynomials evaluated at x = scale as 

(5.27) 7r„(0;t)| t=s/4n ~ (-l) n (a + 3)„ exp f / — C(u)) . 

Proof. The relation (|5.26[) arises from integrating (|3.27p with respect to t and then 
employing the scaling form for the integrand as given by (|5.1ip . The second relation 
is derived by integrating (|3.39|) and using (|5.15p for the scaling of the resulting 
integrand. The factorial and Pochhammer prefactors arise from the normalisations 
at t = 0. □ 

We find that the polynomial ratios have well defined scaling behaviour rather 
than the polynomial themselves. 

Proposition 5.4. The polynomial ratios Q n (x;t), R n (x;t) scale as 

(5.28) lim Qn{x;t)\ x =-z/4n,t=s/4n=<l(z;s), 
n— *oo 

(5.29) lim nR n (x;t)\ x= _ z / 4nt=8 / 4n =p(z;s), 

n — *oo 

where q(z\ s),p(z] s) are entire functions of z. 
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Proof. We start with the product form of the scaled polynomial 



(5.30) 



Qn(~ 



z s , 
An' An' 



AnXj iTl (s/An) 



and seek bounds for the logarithm of the ratio of the n-th scaled polynomial to the 
(n — l)st. When < s, z < oo we can write the general bound as the sum of four 
contributions 



(5.31) 



log 



i(— z/An; s/An) 



Q„_i(-z/4(n-l);«/4(n-l)) 



< 



log 1 



E 

3=1 



log 1 



log ^1 
log ( 1 



\ Anx hn (t) 

z 

A{n-l)x 3 , n (t) 

z 

A(n - l)xj^(t>) 



Anx ntn (t) 
- log f 1 - 



- log 1 



log 1 



4(n - l)a;j,„(i) 
z 

A(n - l)x :)> (t>) 

2; 



A(n - l)xj,„_i(i>) 



where i = s/An and i> = s/A(n — 1). Using the inequality 

>o g i±|<.4- B 
for A > B > we can find a simpler bound 



(5.32) LHS of dOTTl < 



Anx n , n (t) 

n-l 

E 



3=1 



Anxj, n (t) A(n - l)x hn (t) 

z z 



A(n - l)xj, n (t) A(n - l)a; J ,„(i>) 

z 



A(n - l)x 3 . n {t>) A(n - \)xj %n -i{t>) 
Considering the first sum of (|5.32|) we see that this is 



E 

3=1 



Anxj, n (t) A(n - l)xj, n {t) 



An(n 



71—1 

1) ^ j Xi.ni 



< 



< 



7,n (^) 

J = l 



n 1 

1) .. •^■i.n. \ 



4n(n- 1) 2 



zO(n : ), as n — > 00, 
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for all bounded s. The second term of ()5.32|) is 

n-l 

E 



4(n - l)xj>(t) 4(n - l)x,>(i>) 



4(n - 1) ^ 



The summand appearing here can be bounded in the following way 



1 



1 



(5.33) 



< 



I a; j,n(i>) — Xj,n(OI 
%j,n (t> ^)*Ej,n (^) 

-(t>— i) max |ij.„(u)|. 



Now from (|2.96|) we note that 



,. , -1 On(u)+U 
\ X 3,n{u)\ = U — Xj tn (u), 



and so furnishes a bound on (|5 . 33[) 



< 



1 



4d 



max 



?„(u) + u 



< 



Xj,n{t>)xj, n {t) 4n(n — 1) s ue(t,t>) — n (u) + a; J >(w)"' J 

_L 1 , , , On{u) + M 

71 - 1 Xj in (f>)Xj, n (i) ' «£(t,t>) — V n {u) +Xj, n {U) 



%j,n (^) ; 



< 



1 1 



max 



i{u) + u 



This means that 

n-l 



E 

3=1 



1 



n-l Xj>(i>) «e(t,t>) -9 n (u) + x 3 . n {u) 

n (u) + u 



Xj',n(^) Xj. n (t>) 



< 



< 



1 " 
1 \ ~» 

y IT1SX r\ / \ / \ ; 

- 1 «e(t,t>) -0 n (u) + Xj tn {u) x,-,„(t>) 

n-l „ , . n-l 

1 \ - V n (U) + U \ -\ 1 

} max — — y 

- 1 "f— f «e(t,t>) -0 n (u) + Xj.niu) f-j' £j.n(£>) 



3=1 



1 ■r^ 9 n (u) + U 

< > max — — — — y 

n — lj^ «6(t,t>) —OnW + x j,n(u) Xj>( 



< 



< 



< 



E 



max 



(u) +Xj, n (u) ^ a; iin (t>) 

n(«) + W 



9 n (u) + u 



2(n - 1) «e(t,t>) ^ ~6 n (u) + x^ n (u) 



K n (u) - u 



2(n — 1) «e(t,t>) 



?n(u) 



where we have used (|2.94p in the last step. The total contribution of the second 
term is therefore bounded by 



8(n- l) 2 



0(1) = zO^r 1 ). 
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The third sum in ([Q2]) is 

n-1 



4(n- 1) 



2Cj>(i>) Xj,„-i(t>) 



1 



1 



< 



1 



4(n 



E 



~ ^,n-l(*>) 



From the identity (|3.48[) and the scaling of the variables involved as given in (|5 . 1 7[) 
we conclude that the contribution of this sum is bounded by 

-0(1) = zO^" 1 ). 



4(n- 1) 

Finally we note that the isolated term in (|5.32|) is of order zO{n~ 2 ) as a leading 
estimate of the largest zero x n , n is of order 0(n). This establishes that Q n under 
the hard edge scaling of the independent variables converges to a limit as n — > oo, 
for all real, positive and bounded z, s. □ 

The spectral and deformation derivatives of the Q n ,Rn system scale to the 
corresponding derivatives of the q,p system as in the following result. 

Proposition 5.5. Specify scaled quantities as in Proposition I5.il Corollary 1 5. ,21 
and Proposition 15. 21 The spectral derivatives of the q, p system are 



(s — z)zd z q = -zCq — (jit + z)p, 



(5.34) 

(5.35) (s — z)zd z p 

and their deformation derivatives are 



q+ [-2s + z(C + a + 2)}p, 



(5.36) 
(5.37) 



(s — z)sd s q = zCq + (fi + s)p, 

(s - z)sd s p = z£q - [s(2C + a) - zC]p, 



The boundary conditions satisfied by the solutions q(z; s) and p(z; s) of the above 
system on the domain s > z > along z = are 



(5.38) 
(5.39) 

for all s > 0. 



q(0;s) = 1, 
p(0;s) = 0, 



Proof. The first spectral derivative (|5.34[) follows from the scaling of (|3.53p and 
employing (|5. 10|) . (|5.12| and (|5.15j) . The second member (|5.35|) is derived from 
the scaling of (|3.54|) and using (|5.15|) and (|5.17|) . The first deformation derivative 
([536]) follows from the scaling of (pLjg) and utilising ([5~T0)l . (j5T2]) and ([BTTS) . The 
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second deformation derivative (|5.37|) arises from the scaling of (|3.56j) . employing 
()5.16D and (|5.17|) and noting that 



(«n + t)— [K n - t)- 



s(2C + a) 



i a n^oo An 
t=s/4n 

The boundary conditions (|5.38[) and (|5.39|) follow from the definitions (|3.41|) and 
(|3.52|) respectively and the scalings in Proposition 15.41 □ 

Remark 5.1. Both the spectral derivative (|3.43|) and the deformation derivative 
([315]) scale to 

(5.40) (fx + s)zd z q + O + z)sd s q + zC(s)q = 0, 

and this mixed derivative equation can be easily found from Proposition 15.51 by 
eliminating the variable p between (|5.34p and (|5.36|) . i.e. (/i + s) times (|5.34[) plus 
(/i + z) times (|5.36p . The three-term recurrence relation (|3.42p scales to 

(5.41) z 2 d 2 z q + 2szd z d s q + s 2 d 2 s q + s^—^- [zd z q + sd s q] - -^zq = 0. 



This can also be recovered from the equations of Proposition 15. 51 If we eliminate q 
between (|5.34j) and (|5.36|) by adding them we find that 

(5.42) p = zd z q + sd s q. 
Employing this we find that (|5.4ip is equivalent to 

(5.43) zd z p + sd s p + (—1 + s- )p — -zq = 0, 

fX 4 

which can be found by adding (|5.35|) and (|5.37p and noting the relation (|5.20l) . 

Remark 5.2. In addition to the boundary conditions at z = given by (|5.38|) and 
()5.39|) there are also relations along s = z 

d C 

(5.44) —q( S ;s) = —«(«;«), 

as fi + s 

and further 

p(s;s) £ sC 



(5.45) 

q(s; s) C + a Li + s 

However these are a consequence of the spectral and deformation derivatives and 
so do not constitute independent boundary conditions. 

Remark 5.3. The compatibility of the two sets of derivatives, (|5.34p and (|5.35p on 
the one hand, and (|5.36p and (|5.37|) on the other hand, affords a check on our 
results. We find that compatibility of (|5.34p and (|5.36p leads us to conclude that 

(5.46) Z = sC+^(li + s), 
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and we also recover (|5.20j) . Similar considerations applied to ()5.35|) and (|5.37|) imply 
that 

(5.47) s£=-(2C + a + 2)Z+-s(2C + a), 

and we get (|5.46p again. Using (|5.19p to eliminate fi we arrive at a coupled pair of 
first order ordinary differential equations 

(5.48) s c = £ + JLc(C + a), 

(5.49) si = -(2C + a + 2)<e+ ]s(2C + a). 



Using the latter equation to eliminate C we obtain a second order ordinary differ- 
ential equation for £, which by means of the transformation 



(5.50) m = \- SV{3) 



4y(s)-l' 



is transformed into the standard equation for the fifth Painleve transcendent. This 
is a degenerate case of Py which reduces to the third Painleve transcendent be- 
cause the parameters are ay = 9/2, /3y = — a 2 /2,7y = 1/2, 5y = 0. Making an 
independent variable transformation s i— > 2s so that 7y = 1 we determine the 
Pin parameters to be am = 2(2 — a), /3m = 2(a + 3), 7m = 1,5m = — 1 which is 
consistent with those in Proposition 15.21 



The matrix form of the spectral derivatives (|5.34|5.35j) and deformation deriva- 
tives ([5.36I5.37P yield the Lax pairs 

(5.51) d,*={A 00 + ^- + ^-)% 

|_ z — s z J 

(5.52) d s $= is-^-W 



z — s 



iii the matrix variable 



(5.53) tf(z;s) 



g(z;s) 



The system has two regular singularities at z = 0, s and an irregular one at z = 00 
with a Poincare index of 1. This system is essentially equivalent to the isomon- 
odromic system of the fifth Painleve equation but is the degenerate case. The 
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residue matrices are 

(Ms) 



(5.54) A = 

(5.55) As = 

(5.56) Ao = 



s I 

v -2 )' 

f C{s) ^ + - N 
K H(s) ~C( S s)-a, 



Local convergent expansions about the regular singularities take the form for z = 

(5.58) g (z; S )=5>^( S y« + ™, 

m>0 

(5.59) P(*;*)= £ t&(*)* Xo+ro , 

m>Q 

for |z| < s and the initial relations amongst the coefficients are found to be 

(5.60) sxofo = s Xo u o = ~2su°. 

This implies that \o = —2 or xo = and Uq = 0. The latter case applies here as 
both q,p are analytic at z = 0, and in addition we also require p = on z = 0. In 
addition r§ = 1. The recurrence relations for general m are 

(5.61) s(m + 2)u° m = (C + a + m + + Qs - A - ^_ 2 , 

(5.62) smrl = + (m - 1 - - u^. 

For z = s we have the convergent expansion 

(5.63) g (.; S )=^r;(#-^ m , 

m>0 

(5.64) p(z; S )=)]«;( S )( S -^ +m , 

m>0 

for |z — s| < s and where in this case the initial relations are 



(5.65) s(- Xs + Cy = -(// + sK, (C + a + x.)«o=f»^- 

Combining these we get a relation which is identical to the second equality in (|5.45[) 
only if Xs = 0, —a. The former case is the one we must choose as q, p are well-defined 
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and finite on z = s. For general m the recurrence relations are 
(5.66) 

s(C - m)r s m + + s)< = (C - m + l)r^_ x + 
(5.67) 

-CC + s(C + a + m)u s m = - + A + (C + a + m + 1)<_ X + \r a m _ 2 . 

This system has a unique solution only if s 2 m(a + m) 7^ for s > and m > 1 
which in turn means that a ^ -N. We also note that the two sets of coefficients 
are related by 



00 



(5-68) C00 = (-l) m E "~ )s n r° m+n (s), 

and an identical relation for m^(s). 

Proposition 5.6. The determinant D n (x,x) scales as 

(5.69) D n {x,x)\ x= _ z u nt=a / 4!n ~ -4A ri 7r ri (0)7r„ + i(0) [gc^p - p<9 z g] . 

n — >oo 

Proof. We can employ the polynomials Q ra , i? n in the evaluation (|2.45|) and find 

D n (x,x) = A„7r n (0)7r n+ i(0)[Q„ + ii?J i+1 - R n+ iQ' n+1 }. 
Applying the scaling of Proposition 15.41 to this expression we arrive at (|5.69|) . □ 
Proposition 5.7. The distribution A n , a (y) scales to 

(5.70) ^-A n , a (y)\ y=z/4nn ~^A a {z). 

Proof. This is the only possible scaling consistent with the scaling of the indepen- 
dent variable to the hard edge, y = z/4n. □ 

Proposition 5.8. The distribution of the first eigenvalue spacing at the hard edge 
is given by 

z 2 

(5.71) A a (z) = 42Q+3r(a - 1)r(fl - 2)r2(fl - 3) 

ds s a (s - z) a exp ( -- + J —Hv) + 2C(v)] \ [qd z p - pd z q) . 

Proof. We apply the scaling (|5.70p to Eq. (|1.3ip and utilise the previous relation 
for the scaling of the integrand (|5.69|) . For the first three factors on the right-hand 
side of (|5.69|) we can use the scaling results of (|5.26|) and (|5.2T[) , yielding the above 
integral representation. □ 

We note that the factor of the integrand of (|5.7ip can be written as 

5.72) qd z p-pd z q= -q ~-qp-\ p - -{£q - Cp) '—. 

4 z sz E, s — z 
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The corresponding distribution A^(z) defined in (| 1 .3511 and (|1 .36|) for the special 
cases a = ±1/2 is given by 



(5.73) A ± (z) 



4 2a + 2 i> + l)r(a + 2)r 2 (a + 3) 



^ ds s 2a+1 (s - z) 2a+1 (2 S ~ z) 2 exp (~ +1*^-1 



["(«) + 2C(«)] 



x [qd z p~pd z q}\ s _ s 2 



5.2. Special Case a £ Z>o- 

Proposition 5.9. In the special case a £ Z>o we have 
(5.74) 

A n (t)\ t=s/4n ~ £lLti^ ( _l)La/ 2 j (2n) a s -a det[l . 



and 
(5.75) 



(S) = S ^ l0gS adCt ^+ 2 - k ^h.k=l,...,a' 



whilst for a = we find u(s) = 0. 



Proof. The first relation (I5.74[) follows from an application of the Hilb type asymp- 
totic formula (|5.2p to (14. ip and the second follows by using this result in (|3.27p . □ 



Another consequence of the Hilb formula is the scaling of the orthogonal poly- 
nomial ratio as given by (|4.8p . 



Proposition 5.10. The scaled orthogonal polynomial ratio q(z; s) is 



det 



(5.76) q(z; s) = z 



and for a = is 



[(^) fe/2 wv% =1 ,... a+1 

[Ij+2-kiyfs)] , i 

fc=l,...o+l 



det 



[l/8VM/2s,l/s 3/2 ,0,...,0] 

[i>+2-fc(VS)] j=l,...,a 
fe=l,...a+l 



, a > 1 



(5.77) 



<j(z;s) = -I 2 (y/z). 



Finally the eigenvalue spacing distribution (|4. 16[) takes the following form in the 
scaling limit. 
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Proposition 5.11. The distribution of the first eigenvalue spacing at the hard edge 
for a G Z>o is 



(5.78) A a {z) = 2-^J ^ 



e -./4 



x det 



[Ij+2-k(Vs)] j=l,...,a 

[(s/z)^- k y^ j+2 . k (^)} +2 j=1 . 



fe=l,...a+2 



for a > 1 and for a = is 

(5.79) A Q (z) = ± e -*/ 4 - /x(^)7 3 (^)] . 

Proo/. We apply the Hilb asymptotic formula (JOJ) to (|3~T5|l . As |(a + l)(a + 2) + 
1 + [a/2j is always even for a £ Z we have (|5.78|) . □ 

6. Analytical Studies at the Hard Edge 

In this section of our study we intend to develop the analytical and non-formal 
theory of the solutions to the defining ordinary and partial differential equations 
described in the previous sections. This is because we wish to compute precision 
numerical data characterising the distribution function of the first eigenvalue spac- 
ing at the hard edge A a (z) for arbitrary parameter a. For this purpose it is not 
sufficient to employ a single local expansion of the cr-function, about s = say, 
because it has a finite convergence domain and one cannot use this to evaluate the 
s-integrals on the interval [0, oo). For this reason we construct a patchwork of over- 
lapping local expansions including Taylor series expansions about regular points 
So for positive and real values. A similar approach was undertaken by Prahofcr 
and Spohn in their study [29] of the exact scaling functions for one-dimensional 
stationary KPZ growth. 

6.1. The cr-function expansion about s = 0. The terms given in the expansion 
of the Painlcvc III' cr-function (|5.22p are the minimum required to specify the full 
non-analytic Puiseux-type expansion of the particular solution for v(s) about s = 
in the sector —it < arg(s) < ir. This is the primary data specifying our particular 
solution and we need to use this in various ways in order to compute the distribution 
A a {z). 

Proposition 6.1 (|18j.[2Tj). The Painleve III' a -function u(s) has a Puiseux-type 
expansion about the fixed regular singular point s = of the form 



oo oo 

x x „j+ka 



(6.1) v(s) = 2^ 2^ c fcj s J+fca , a G C, < Re(a) < 1, 

fc=0 j=0 

with |arg(,s)| < ir and is convergent in a finite domain s G {z g C : \z\ < R, \z a \ < 

Ra}- 
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The coefficients Ckj are determined by recurrences which follow from the sub- 
stitution of expansion (|6.1|) into the relation (|5.25|) . We will assume that a is not 
rational for simplicity. Considering terms in the resulting equation with s p and 
p G Z>o then the p = case implies that if cq.q = then this choice fixes 

< 6 - 2 ) ^ = At^y 

The p = 1 is automatically zero but for the p = 2 case one has the two options 

(6 ' 3) C °^ 8(a + 3)( a + 2) 2 (a + ir " °' 

The former case applies here by comparison with (|5.22p . For p > 3 the recurrence 

is 

p— l 

l a(p + a+l)(p-fl-3) 2 -i '\ 

(6 ' 4) ~ 2 (a + 3)(a + 2)2(a + l) ^ " ^T2 ^ j(P " 

p-i 

+ 8c ,2 2j(j - 1)(P - j)co,j-lCo,p + l-j 

P-l 

+ 53j(p + 2 - i)[(j - l)(p + 1 - j) - (a + 2) 2 ]c 0j co,p + 2- J 

p-lp+2-j 

+ 4^ ^ 3 m (P+ 1 - J ~ "l)cojCo,,„C 0i p + 2-j-m = 0, 
j=3 m=l 

which allows for co lP to be recursively found, valid for a £ Z. Terms with s ^ 2 +9 a for 
q = 2k > 2 imply that c^o = given that co,o = and a^3/ (fc — 1), — l/(fc + 1). 
Consequently all terms with s ~ 1+<?a vanish also. The generic recurrence relation 
following from examining the s p+qa term is 

(6.5) 

i 9 P+l 

-a(a + 2)(p + 1 + qa)c qtP+1 = ^(j + fca)[p - j + (g - k)a]c k ,jC g - k ,p+i-j 

k=0 j=0 

q p+2 

-^^(i+ka)[p+2-j+{q-k)a][{j-l+ka)[p+l-j+(q-k)a^ 

k=0 j=0 

q q—kp+2p+2—j 

X] (.7+ fca )( m +^ a )b+ 1 -J- TO +(9-fc-0 a ]CfcjQ,mCg_fc-i.p+2-j-m, 
k=0 1=0 j=0 rn=0 

for q,p > 0. The convention is taken that sums with upper limits less than their 
lower limits are zero, or equivalently coefficients with negative indices are zero. 
Contrary to appearances the highest coefficient in (|6.5[) will turn out to be c q%p as 
Cq,p+2 occurs with co,o as a factor and c 9iP +i has a factor of 

(6.6) -a(a + 2) + c , - 2(a + 2) 2 c oa - 8c 0i0 c ,i, 
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which is zero as a consequence of co,o = and (|6.2p . We find that the coefficient 



Cq :P occurs as a linear term and has a factor of 

(6.7) (p-l+ga)c ,i(4c ,i-l)+4(]3+ga)(p-l+ga-(a+2) 2 )co,2-16(p+9a)co,oCo,2- 

With the above evaluations (|6.2p and (|6.3[) this is 

, v _ o(p + 1 + (g + l)o)(p - 3 + (g - l)a) 

1 ' ' 2(a + 3)(a + 2) 2 (a + l) 

This vanishes when q — l,p — 3 and wc find that C1.3 is undetermined and therefore 
is a free parameter. With these choices a triangular subset of the coefficients vanish 

(6.9) Cfcj =0, j = 0, . . . , 3k — 1 for k = 1, . . . , 00, 

so that the initial non-zero term in the j— sum is cu.zk- All other terms are fixed 
by ci,3 and a. 

For the variable /x(s) we can deduce the following Puiseux-type expansion from 
the above work 

(6 ' 10) ^ S) sZ -^h S - (a + 3)(a + 2na + l/ 

2 -s a + 3 . 



4 a + 2 (a + 2)(a+ l)r 2 (a + 3) 

In regard to the quantity 2C(s) we require as much detail about this as for v(s). 
This variable is a cr-function for an identical problem, where the fixed exponent 2 
in the weight (|1.32p is replaced by 3. If we define the expansion coefficients for this 
object 

00 00 

(6.11) 2 C( s ) = EE at .J s3+ ' :a ' aGC ' 0<Re(a)<l, 

k=0 j=0 

then the defining recurrences for these using (|5.20j) arc 

(6.12) a , = 0, a ,i =4(a + 2)(a+l)c , 2 , 
and for j > 2 



(6.13) a ,j = 2(a + 2) - (j + l)(j - 2 - a)c j +1 + lapj+i- 



lCo,l 



i=i 



For the case k > 1 we have a^,j = for j = 0, .., 3A; — 1 and the remaining non-zero 
terms are given by 



(6.14) oij = 2(o + 2) - (j + l + a)(j -2)ci 



3+1 



j 3+1 
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for k = 1 and the general case k > 2 by 
(6.15) a k ,j = 2(a + 2) [ - (j + 1 + ka)(j -2 + (k- l)a)c K 



i=0 



+ fca ) a 0,i + l-/Cfc,i + ^ £afc,j+l-ZQU 
Z=2 

fc-i i+i 

+ E E^ + rna)ak- m ,j+i-iCm,i 



m=l 1=0 



The first few terms are 
(6.16) 2C(s 



a{a 2 - 5a - 18) 



s~*o 2(a + 3)(a + 2) 



+ 



8(a + 4)(a + 3) 2 (a + 2) 2 (a+ 1) 
1 

4 a + 2 (a + 2) (a + 1)I> + 3)r(a + 4)' 



„a+3 



and 

(6.17) £( s ) 



3o 



^o4(a+3) 8(a + 4)(a + 3) 2 (a + 2) 



s a+4 



4 a + 7 (a + 3) (a + 2) (a + l)r 2 (a + 4) 

6.2. The g,p expansion about s = 0. In this subsection we seek local expansions 
about s = for the coefficients functions ^(3)^^(3), r^^^u^s) appearing in 
(|5.58p . (|5.59|) . (|5.63p and (|5.64|) . In parallel with the transcendent quantities these 
will have Puisuex-type expansions of the form 



(6.18) 
(6.19) 



r m( s ) — E r m,k > j S ^ +ka \ 
k.j>0 

J+ka 



.(-) = E 



fc.j>0 



with analogous expansions for the remaining two coefficients. This is immediately 
clear because the recurrence relations for these imply that they are polynomial 
functions of the variables n(s),C(s)^(s). These recurrence relations imply the 
following ones for the z = coefficients r® m k ^ , k j 

1 



(a + 171 + 2)u rn ^ kQ — -^ r m-l,k,0> mr m.k,0 — u m,k,0' 



(6.20) 

for m > 1, k > 0. For the general case j > 1, k > we have 

^ k j 

k 3 



^ k j 

(6.21) (m + 2)u° mXj _ 1 = (a + m + l)u° m _ hktj + ? E E 



a k-q,j-p 



9 =0 p=0 



5=0 p=0 



1 



C k-q,j-p + r, a k-q,j-p 





m— l,Q,p 

r o _1 

m— 1 ,q,p ^ m — 2,k,j ' 
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and 



(6.22) mrl^, = (m-l)r° m _ 1Aj ~ \ £ £ - 



a k-q,j-p r m -l,q,p 



2 

9 =o p=0 



fc i 

q,j—p u rn,q,p ~ U m—X,k,jl 



+ u m,k,j-i - 4 5Z ~~ p + ( fc ~ g)a]c fe _„.,_„ui „ „ - U° 



g=0 p=0 

These equations can be solved for successive values of m starting with the m = 
values for all k,j 

(6.23) u o,k,j = 0' r o,k.j = 0, except for k = j = where rj] = 1> 

which follow from r§(s) = l,u{j(s) = 0. The next few coefficients can be read off 
from 

(6.24) r?( S ) 



s -o4(a + 3) 16(a + 4)(a + 3) 2 (a + 2) 

1 



128.4 a (a + 3) (a + 2)(a + l)r(a + 5)r(a + 4) 



(6.25) u?(«) 



s -o4(a + 3) 8(a + 4)(a + 3) 2 (a + 2)' 

1 



128.4 Q (a + 3) (a + 2) (a + l)r 2 (a + 4) 



s a + 3 , 



(6.26) r§(«) 



32(a + 4)(a + 3) 64(a + 5) (a + 4) (a + 3) 2 (a + 2) 

1 



512.4 a (a + 3) (a + 2) (a + 1)I> + 6)T(a + 4) 



s a+ \ 



(6.27) u° 2 (s) 



-*o 16(a + 4) (a + 3) 64(a + 5)(a + 4) (a + 3) 2 (a + 2) 

1 



512.4 a (a + 3) (a + 2) (a + l)r(a + 5)r(a + 4) 
The analogous results for fc . , u s m k ■ are 

(6.28) (a + m + 2X, jfc)0 = -^r s m _ lkfi , mr s mkfi = u^ A0 , 



s a + 3 . 
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for m > 1, k > 0. Again for the general case j > 1, k > we have 

1 fe j'-i 

(6.29) (a + mX^j-i + - ^ ^ afe- 9 J-i-p«m, 9 , P 



g=0 p=0 



^ k j 

(a + m + l)u s m _ ljkjj + - ^ Yl 



r> / , a k-q,j-p U m—~L,q,p 
" q=0 p=0 



1 



k j-1 

q=0 p=0 

k J 

£ 53? -J>+ (*-«)<»] 

q=0 p=0 



c k-q,j-l-p + e^k-q^-l-p 



m— l,fe,j — 1 



C k-q,j-p + n a k-q,j-p 



m—l,q,p 



m-2,k : ji 



and 



fc j-l 



(6.30) mr^j^ - - ^ ^ afc- g ,j-i 



9 =0 p=0 



fe i 



2 /L~i ak -1^-p r m~l,q,p 

fe j 



U m-l,fcj + 4 5Z ~ P + ( fc ~ ^alCk-qj-pU^ p, 

g=0 p=0 

Again these can be solved for successive values of m starting with the initial values 
given by the relations 



(6.31) 



'0,k,j 



n=0 



along with an identical formula for Uq k ■ . 

6.3. The (T-function expansion about ,s = oo. The nature of the expansion of 
v(s) about s = oo is rather different because this fixed singular point is irregular 
in the case of the third Painlcve transcendent. 

Proposition 6.2. The formal asymptotic expansion of v(s) about s = oo has the 
form 



(6.32) 



where d-\ = ±^a. 



-j/2 



j'=-i 



Proof. We start with the general ansatz of 

(6.33) v{s) =als ka +0{s {k - 1)a ), 
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where k £ N and a £ C with < Re(a) < 1. Using (|5.25j) we find that the only 
terms which can balance the 0(1) term are the 0(s 2fca_1 ) and 0(s 3ka ~ 2 ) terms. If 
ka < 1 then the former choice applies and we find ka = 1/2. If we assume ka > 1 
then the latter case must be chosen but this leads to ka = 2/3 in contradiction to 
the hypothesis. With the correct choice of ka we find the above relation for leading 
coefficient, Considering the sub- leading terms we note that it is only possible 
for the terms of 0(s~ Q ) to balance those of 0(s -1 / 2 ) so that in fact a = 1/2 
and k = 1. This is entirely consistent with the expansion for the Painleve III 
transcendent qm(t) about t = oo as found in |15j . □ 

The choice of the sign of the leading coefficient is positive in our application. 
Consequently we find the first few terms of the asymptotic expansions of the a- 
function 

(6.34) V {a) = loa^-iaCa + ^ + H^-i/a + HaV 1 

s^oo 2 4 lb lb 

+ ^a(16a 2 -7) S - 3 / 2 + 0( S - 2 ), 
25b 

and the auxiliary variables 

(6.35) fi(s) + s = as 1 / 2 - H as -V2 _ 15 a 2 s -i 



S— OG 



4 



^U(16a 2 -7) S - 3 / 2 + 0( S - 2 ), 



and 



(6.36) 2C(s)+a = las- 1 ' 2 + + ^-a(8a 2 - 21)s~ 3 / 2 + 0(«- 2 ), 

s^oo 2 2 16 

and 

(6.37) e( S ) g = x i«V a - go," 1 " - fa 2 ,- 1 + 0(s-^). 

The large s-regime also implies a simplification of spectral derivative (|5.51[) 
which becomes 

a 1 



(6.38) d z *=[ 2 ^- s ) a Z 2 {l + 0( S -V^)} 

Using the substitution \I> i— > exp(— a ln(s — z)l/2)^l this decouples and can be solved 
in terms of the modified Besscl functions. An application of the boundary condition 
(|5.38jl implies that 

(6.39) q (z- 8 )„(-LS\ 1 *I 2 {Vz), 

\s — z J z 

(6.40) p(z;s) ~ 7 z-f f-J 3 (V5) 

V s — z J az V z 
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6.4. The c-function expansion about a regular point. In this subsection we 
seek Taylor series expansions for the sigma function and derived variables about 
regular points So, taken to be positive and real without any loss of generality. Let 
us write 

oo 

(6.41) v(s)=J2<ij(s-so) j , 

3=0 

and using (|5.25j) we find the following recurrence relations for the coefficients dj 

(6.42) Asldj - ^a 2 + -a(a + 2)dt + d d l - [s + (a + 2) 2 ]d\ + 4(s di - d )d\ = 0, 
and for n > 1, d 2 ^ 

(6.43) 4s 2 (n + 2)(n + l)d 2 d n+2 = ~a(a + 2)(n + l)d n+1 

n 

+ [s + (a + 2) 2 ] Y,U + - J + l)d j+1 d n - j+ i 

3=0 

n 

+ + - j - l)rfj+id„_ 3 
3=0 

re— 1 

- + !)i( n - i + ~ j)dj+id„- j+ i 

3=1 

n-1 

- 2s ^(j + 2)(j + l)(n - j + l)(n - ])d J+2 d n _ J+1 

3=0 

n-1 

- s o E^' + 2)0" + l)(n -i + 2)(n - j + l)d J+2 d„_ J+2 



3 = 1 
n n—i 



- 4^ 2j(i + l)(j + l)(n - i - j - l)d i+ idj+id„_i_j 

n re— i 

- 4s ^^(i + l)(j + l)(n - i -j + l)d i+1 dj + id n -i- j+ i. 

i=0 j'=0 

These recurrences are solved subject to the initial values of 

oo oo 

(6-44) do = EE c ^ s 5 +fc °> 

fe=0 j=3fe 
oo oo 

(6.45) d i = E E + Mc fc ,3-4" 1+fa > 

fc=0 j=3fc 



PETER J. FORRESTER AND NICHOLAS S. WITTE 



which in turn can be found from the solutions to the recurrences (16.41) and (16.5 



In addition we define 








(6.46) 


fx(s) 


i>o 




(6.47) 


C(s) 


j>0 


soV 


(6.48) 




3>0 





The coefficients appearing here are computed using the recurrences 
(6.49) 

/o = «o(4di - 1), fi = 4di - 1 + 8s d 2 , /i=4[ 7 'd i + so(7 + lK-+i]. J > 2, 
and subject to /o ^ 

(6.50) 2/o<?o - -2s - (a + 3)/ + s /i, 

(6.51) 2/o. 9 i = -2 - (a + 2)/j + 2s / 2 - 2/i<7o, 

i-i 

(6.52) 2/oft = (j -a- 3)/j + (j + l)a f j+1 - 2 £ /,„ fcfffc , j > 2, 

and provided so + /o 7^ 

(6.53) (s + f )h = -s Q g (a + g ), 

(6.54) (s + fo)hi = -as gi - 2s .9o3i - (a + ffo).9o - (1 + /i)^o, 

(6.55) (s + f )hj = -as gj - ag J ^ 1 - s ^ gj-k9k - £ 9j-k9k-i 

k=0 k=l 

i-2 

-(I + WVi-^/h^, i>2. 

fc=0 

6.5. The expansion about a regular point. We will seek a Taylor series ap- 
proximation for the scaled polynomial and associated function q(z; s),p(z; s) about 
the regular point (zq, sq) with < zq < sq. Let us write 

(6.56) q(z; s) = J2 r i,fc(-2 - ^) 3 {s - s Q )\ 

j,k>0 

(6.57) p(z;s)= ^2 u jyk (z - z ) J (s - s ) k , 

j,k>0 

Using the first spectral derivative (|5.34j) we obtain the recurrence relation 
(6.58) 

zo(so-zo){j + l)r j+ i ! k + (so-2zo)jr jt k-(j -1^3-1^ + z (j +1)^+1^-1+^^-1 

k 

= -zou j:k - Uj-i t k - £ [fk-iUj t i + gk-irj-\.i + z gk-irj t i] , 

2=0 
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for j,k>l. When k = and j > 1 we have the specialisation 

(6.59) z (s - z )(j + 1)^+1,0 + [(so - 2z )j + z go\r jt o + [g - j + lfo-1,0 

= -(/o + zo)uj,o ~ Wj-i,o- 

The second spectral derivative (|5 . 35[) yields the recurrence relation 
(6.60) 

z (s - z )(j + l)uj + \ t k + [(s - 2z ).7 + 2s - (a + 2)z ]u Jtk - (a + j + l)uj-_i,fc 

+ zo(j + l)«j+i,fc-i + (j + 2)Uj,k-i 

1 . , 1, , 1 1 1 

= ^o(s - «o)rj,fc + -(s - 2zo)rj-i,k - -^ r j-2,k + -zor^k-i + -rj-i, k -i 

k k 

- ^2 ft.fc-i( r j-i,i + z o r j,i) + y^9k-l(uj-i,i + z Uj.i), 

1=0 1=0 

for j, k>l. When k = and j > 1 we have the specialisation 
(6.61) 

zo(so-zo)(j+l)uj+ifi+[(so-2zo)j+2so—(a+2)zo-zogo]uj,o—[a+l+go+j]uj-i t o 
= [^-zo(so - zo) - hozolrjfi + [^(s - 2z ) - ft-o]^-i,o - ^ r j-2,o- 
Using the first deformation derivative (|5.36p we obtain the recurrence relation 
(6.62) 

so(so-zo)(k+l)rj t k+i-so(k+l)rj-i t k+i + (2so-zo)krj tk -krj-i ! k+(k-l)rj t k-i 

k 

= sou J:k + uj.k-i + ^ [fk-n,i + gk-irj-\,i + z gk-irj,i] , 

1=0 

for j, k > 1. When j = and k > 1 we have the specialisation 

(6.63) s (s - z )(k + l)r ,fc+i + (2s - z )kr ,k + (k- l)r ,k-i 

k 

= s uo,fc + wo,fc-i + y~"}fk-iUo,t + Zo9k-iro,i\- 

1=0 

The second deformation derivative (|5.37l) in turn gives us the recurrence relation 

(6.64) s (s - z )(k + l)uj,fc+i - s (fc + l)u 3 -_i,fc + i 

+ [(2s - z )k + as ]uj } k - fritj-i.k + (a + fc - l)uj,fc_i 

A! A? A.* X 

= ^ ftfc-j( r j-l,J+ 2; Or3,i) _ ( 2s O-«o) gfc-i^j.i+y^ gk-lUj-l,l-2 ^ 9k-l-lUj,l, 
1=0 1=0 1=0 1=0 
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for j, k > 1. When j = and k > 1 we have the specialisation 

(6.65) s (s - z ){k + l)u .fe+i + [(2s - z )fc + aso]Mo,fc + [a + fc - l]«o,fc-i 

fc k fc-1 

= z ^ h k _ir 0j i - (2s - z ) ^ 9k-lUo,i - 2 ^ 9k-i-iu 0: i- 

1=0 1=0 1=0 

These recurrences can be solved in the following way. First (|6.59[) and (|6.6ip arc 
solved for TjQ,Ujo f° r J > 1 in terms of roo,Uoo. Then these solutions can be 
substituted into the boundary conditions 

(6.66) XV j)0 (-«>) J ' = l, 

j>0 

(6.67) $>j,o(-*o) J = 0, 

i>o 

and this allows for tq.o, ito,o to be found. Next the sequence ro.fc, ito,fc can be found 
for k > 1 using (|6.63| and (|6.65p . Finally the two general systems (|6.58p . (|6.60p and 
(|6.62p . (|6.64[) can be employed to compute r^k^Uj^k for j, k > 1. 
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6.6. Numerical studies at the Hard Edge. Using the integral formula for the 
distribution A a (z) as given by (|5.7ip it is possible to compute values of this and 
the examples of a = 0, 1, 2 are plotted in Figure 2. However we wish to characterise 
it in a precise quantitative way and evaluate the moments of this distribution 

/>oo 

(6.68) m k := / dz z k A a (z), k G Z> . 

Jo 

These are easily seen to be 

(6.69) m k = 42a+3r(a - 1)r(fl - 2)r2(fl - 

pOO pi 

X / ds s k+2a+3 e -s/4 e F(s) / ^ ^+2 (l _ ^aQ^. ^ 

where 

(6.70) F(s):= f— [u{v)+2C{v)l 

Jo v 

and 

(6.71) G(z;s):=qd z p-pd z q. 

We note that by employing the large s asymptotic form of q(z; s),p(z; s) as given in 
(|6.39I6.40[) we can deduce the asymptotic form of the spacing distribution is given 

by 

(6.72) A a ( Z ) ~ Ce -z/±+(a+2)^ z -l/2-*V\ 

z — 'OO 

where C is a constant which cannot be found from our methods. For small z we 
find that 

(6.73) A a (z) 42a+3r(a - 1)r(a - 2)r2(fl - 3) 

*» 2 °(b-^)^(-MvM») + ^(»)i)« 2 . 

where we have used (|5.58I5.59[) and the fact that uj(s) = 1/12 — ^(s)/3s. Therefore 
we can conclude that the moments exist for Re(/c) > —3, Re(a) > —1 and Ke(k + 
2a) > —4. An instance where exact evaluation of the moments can be made is the 
case a = and the first four of these are 

(6.74) m! =4e 2 [/ (2)-7 1 (2)], m 2 = 32e 2 / (2), 

(6.75) m 3 = 384e 2 [2/ (2)+/i(2)], m 4 = 2048e 2 [13/ (2) + 9Ji(2)] . 

We investigated the distributions A a (z) and A ± (z) for the two special cases of 
a = ±1/2 in some detail because of the motivations provided by (|1.6[) and p.7[) . 
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The analogue of (|6~69"1) for A ± (z) is 

1 



(6.76) m 



± 



2 4a + 5 r(a + l)r(a + 2)r 2 (a + 3) 
l 



ds s k/2+2a+3 e - s /4 e F(s) / du u ^ {i - {2 - u) 2 G{u{2 - u)s; s), 



for a = ±i. The statistical data for A a (z) for the cases a = —1/2, 0, 1/2, 1, 2 are 
given in Table 1 and the data for A ± (z) is given in Table 2. 

Our strategy is that by employing local Puiseux-type and Taylor expansions for 
the two factors in the integrand, namely e F ^ and G(z;s), within a given finite 
member of the patchwork of local expansion domains the above integrals restricted 
to this domain can be exactly evaluated. This is essential as numerical quadrature 
algorithms implemented in cither computer algebra software or compiled language 
packages (e.g. QUADPACK) have minimal attained error tolerances which cannot 
be reduced below a fixed bound. For the compiled language option with a floating 
point representation of 64 bits the best one could expect is a relative error of around 
10~ 15 but often it is far worse and around 10~ 8 — 10~ 9 . To illustrate this we have 
computed the statistical data for the a = 1 , 2 cases using QUADPACK routines 
and the results are displayed in Table 1 . In the case of the Puiseux-type expansions 
the integrals are 

(6.77) [ S ds s fe+3+i+n+(m+2)a e - S /4 ) 

JO 

and 

,1/2 ,1 

(6.78) / duu k+2+n (l-u) a or / du u k+2 {l - u) a+n , 

JO Jl/2 

which for k,l,m,n £ Z>o can be evaluated in terms of radicals, the Gamma func- 
tion at integer arguments, a rational function of a and the Whittaker function 
M(a, /3; 5/4) or its specialisations depending on a. For example in the case a = 1/2 
the s-integral reduces to the error and exponential functions. For the Taylor ex- 
pansion case we have the double integral 

(6.79) f 2 ds s k+3+2a (s - s ) m e~ s/4 I du u k+2 (l - u) a (us - z ) 1 , 

Js! Jo 

which for k,l,m € Z>o and so € (si, S2) is evaluated in terms of a rational function 
of a, a polynomial function of sq,z and the Whittaker functions M(a, f3; Si^/4). 
Therefore the only sources of error are from the truncation of the expansions and the 
finite number of intervals, both of which can be adjusted to reduce the contributing 
errors. 

The computations were performed, in most cases, using the computer algebra 
system Maple with a sufficiently large number of decimal digits and found that 
250 digits was more than adequate. In addition we found that we had to tailor 
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the numerical parameters for each case of a = ±1/2 differently as the errors varied 
quite strongly with a (this was especially pronounced as a approached —1). We 
discuss the case of a = 1/2 first. In regard to the truncations about the singular 
point s = 0, we found that 1890 terms in the expansion of the transcendent vari- 
able (|6.1| with k < 30, j + k < 120 yielded an error for the second derivative of 
v{s) at s — 2 which was estimated to be 4.6 x 10 -98 . At most 100 terms were re- 
tained in each of the expansions of the transcendent variables about regular points 
(f6T4T|) . (fOS"|) . (frT47|) and (f6T48|) because much fewer, of the order of j, k < 20, were 
required in the corresponding expansions of the linear variables. The number of 
intervals in the s-direction was taken to be 19 with the sequence of so values being 
{0, 2,5/2, 3, 4, 6, 9, 13, 19, 25, 30, 38, 54, 72, 90, 115, 150, 200, 300, 500}. The bound- 
aries of the s-intcrval with node sq were taken to be located at the midpoints 
of so arL d its neighbouring nodes. This sequence of nodes was chosen to be close to 
an optimal situation yielding the largest separation of each node from its preceding 
node, yet close enough so to ensure that the error in v"{s) at the node was less 
than 9.8 x 10 -37 . For each s-interval with node so two expansion points in the 
z-direction were chosen because a single expansion point could never ensure that 
all of the integration region would fall within the domain of convergence about 
that point. The two points that together yielded the largest convergence domain 
were found to be located at z$ = and zq = sq. Subdividing the z-interval into 
three sub-intervals was found to contribute a variation of less than 3.2 x 10~ 19 
to the normalisation. Another criteria that the sequence of sq nodes had to sat- 
isfy was that each (s,z) integration region fell completely within the union of the 
convergence domains about (so,0) and (sq,sq). For the expansions of the linear 
variables about the lines z = 0, s and about the singular point s = 0, as defined 
in (|6.18I6.19P along with (|5.58I5.59I5.63I5.64|) . we chose the cut-off in the sum to 
be 20. The expansions of linear variables about the lines z = 0, s and about a 
regular point s = so > as defined in (|6.56|6.57j) were cutoff at 25. An overall es- 
timate of the accuracy is provided by the normalisation, which was unity to within 
1.6 x 10~ 18 . The second case with a = — 1/2 was more demanding computationally. 
We needed 5150 terms in the expansion of the transcendent variable (|6.ip with 
k < 50, j + k < 200 and computed these with compiled code using the multiple- 
precision library MPFUN [3], [2]. The error for v" '(s) at s = 2 was estimated to be 
around 7.9 x 10" 59 . Again only 100 terms were retained in each of the expansions of 
the transcendent variables about regular points. A larger number of intervals in the 
s-direction were employed, namely 24, and the sequence of so values was taken to be 
{0, 2,5/2, 3, 7/2, 4, 5, 6, 7, 9, 11, 14, 18, 22, 28, 36, 46, 58, 72, 90, 114, 144, 180, 220, 300}. 
This time the sequence of nodes was chosen to ensure that the error in v"(s) at 
each node was less than 3.6 x 10 -60 . And again for each s-interval with node sq two 
expansion points in the z-direction were chosen at zq = 0, sq. In the expansions of 
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the linear variables about the lines z = 0, s and about the singular point s = the 
cut-off in the sum was chosen to be 20 as before. The expansions of these variables 
about the lines z = 0, s and about a regular point s — sq > was terminated at 
the cutoff of 25 also. The estimate of the accuracy provided by the normalisation 
was 4.6 x KT 20 . 

Because the raw moments grow rapidly with order we have computed some 
standard statistical quantities instead using the definitions of the variance c 2 , the 
skewness 71 and the kurtosis excess 72 

(6.80) a = n 2 , 71 = -— , 72 = — -3, 

in terms of the central moments 

(6.81) fi 2 = m 2 - ml, 

(6.82) ^3 = 7713 — 3toiTO2 + 2mj, 

(6.83) /i4 = ttia — + Qm\m2 — im\. 
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Figure 2 . The distribution of the first eigenvalue spacing at the 
hard edge of random hermitian matrices A a (z) for integral values 
of a = 0,1, 2. 
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Table 1. Low order statistics of the distribution A a (z) for various 
values of the parameter a. 
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